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Abstract 

We investigate the asymptotic zero distribution of Heine-Stieltjes polynomials - polynomial 
solutions of a second order differential equations with complex polynomial coefficients. In the 
case when all zeros of the leading coefficients are all real, zeros of the Heine-Stieltjes polynomials 
were interpreted by Stieltjes as discrete distributions minimizing an energy functional. In a 
general complex situation one deals instead with a critical point of the energy. We introduce the 
notion of discrete and continuous critical measures (saddle points of the weighted logarithmic 
energy on the plane), and prove that a weak-* limit of a sequence of discrete critical measures is a 
continuous critical measure. Thus, the limit zero distributions of the Heine-Stieltjes polynomials 
are given by continuous critical measures. We give a detailed description of such measures, 
showing their connections with quadratic differentials. In doing that, we obtain some results 
on the global structure of rational quadratic differentials on the Riemann sphere that have an 
independent interest. 

The problem has a rich variety of connections with other fields of analysis, some of them are 
briefly mentioned in the paper. 
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1 Generalized Lame equation 

Let us start with a classical problem more than 125 year old. Given a set of pairwise distinct points 
fixed on the complex plane C, 

A = {ao,ai, . . . ,ap}, (1.1) 

(p G N), and two polynomials, 
p 

A{z) = Y\iz - ai) , ^(z) = az^ + lower degree terms G Pp , a G C, (1.2) 

i=0 

where we denote by P„ the set of all algebraic polynomials of degree < n, we are interested in the 
polynomial solutions of the generalized Lame differential equation (in algebraic form), 

A{z) y"{z) + B{z) y\z) -n{n + a- l)K(z) y{z) = 0, (1.3) 

where Vn is a polynomial (in general, depending on n) of degree < p — 1; if deg V = p — 1, then V is 
monic. An alternative perspective on the same problem can be stated in terms of the second order 
differential operator 

C[y]{z)':^A{z)y"{z) + B{z)y'{z), 
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and the associated generalized spectral problem (or multiparameter eigenvalue problem, see [80]), 



C[y]{z) = n{n + a - l)Vniz) y{z) , n G N, 



(1.4) 



where Vn G Pp-i is the "spectral polynomial". 



Special instances of equation (1.3) are well known. For instance, p = 1 corresponds to the 
hypergeometric differential equation. Case p = 2 was studied by Lame in the 1830s in the special 
setting B = A' /2 , aj £ M, and oq + ai + 03 = 0, in connection with the separation of variables in 
the Laplace equation using elliptical coordinates (see e.g. [Ml Ch. 23]). For the general situation of 
p = 2 we get the Heun's equation, which still attracts interest and poses open questions (see [60]). 



Recently, equation ( 1.3 ) has also found other applications in studies as diverse as the construction 



of ellipsoidal and sphero-conal /i-harmonics of the DunklLaplacian [811182]. the quantum asymmetric 
top [H |9l |25], or certain quantum completely integrable system called the generalized Gaudin spin 
chains [23], and their thermodynamic limits. 

Heine [29] proved that for every n G N there exist at most 



a{n) 



n + p 
n 
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(1.5) 



different polynomials Vn such that (1.3) (or (1.4)) admits a polynomial solution y = Qn G IPn- These 
particular Vn are called Van Vleck polynomials, and the corresponding polynomial solutions y = Qn 
are known as Heine- Stieltjes (or simply Stieltjes) polynomials. 

Heine's theorem states that if the polynomials A and B are algebraically independent (that is, 
they do not satisfy any algebraic equation with integer coefficients) then for any n G N there exist 
exactly (j{n) Van Vleck polynomials Vn, their degree is exactly p — 1, and for each Vn equation 
(1.3) has a unique (up to a constant factor) solution y of degree n. The condition of algebraic 
independence of A and B is sufficient but not necessary. It should be noted that the original 
argument of Heine is far from clear, and even Szego [75] cites his result in a rather ambiguous 
form. Recently a significant research on the algebraic aspects of this theory has been carried out 
by B. Shapiro in ^66j, and we refer the reader to his work for further details. In particular, it has 
been proved in [66] that for any polynomials A and B like in ( 1.2 ) there exists G N such that for 



any n > N, there exist fT(n) Van Vleck polynomials Vn of degree exactly p 
a polynomial solution of degree exactly n. 



1 such that (1.3) has 



Stieltjes discovered an electrostatic interpretation of zeros of the polynomials discussed in 

;ntion to the problem. He 
and that all residues pk in 



which attracted common attention to the problem. He studied the problem (1.3) in a particular 
setting, assuming that A C 



B{x) 



k=0 



Pk 



(1.6) 



are strictly positive (which is equivalent to the assumption that the zeros of A alternate with those 
of B and that the leading coefficient of B is positive). He proved in |73j (see also |75l Theorem 6.8]) 
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that in this case for each n G N there are exactly cr(n) different Van Vleck polynomials of degree 
p — 1 and the same number of corresponding Heine-Stieltjes polynomials y of degree n, given by 
all possible ways how the n zeros of y can be distributed in the p open intervals defined by A (see 
Section [2]). Further generalizations of the work of Heine and Stieltjes followed several paths; we will 
mention only some of them. First, under Stieltjes' assumptions C M and pk > 0), Van Vleck |78j 
and Bocher jH] proved that the zeros of each Vn belong to the convex hull of A (see also the work 
of Shah [ 62| [631 EH ES] ) . Polya [5T showed that this is true for ^ C C if we keep the assumption of 
positivity of the residues /O^. Marden [42j, and later, Al-Rashed, Alam and Zaheer (see [2]-[3], |84j . 
[85]) established further results on location of the zeros of the Heine-Stieltjes polynomials under 
weaker conditions on the coefficients A and B of (1.3). An electrostatic interpretation of these zeros 
in cases when ^ C M and some residues pk are negative has been studied by Griinbaum [27], and 
by Dimitrov and Van Assche [16]. For some interlacing properties, see e.g. jlOj . 

We are interested in the asymptotic regime (so called semiclassical asymptotics) when n (the 
degree of the Heine-Stieltjes polynomials) tends to infinity. First general result in this direction, 
based precisely on the Stieltjes model, is due to Martmez-Finkelshtein and Saff [l5]. There the limit 
distribution of zeros of Heine-Stieltjes polynomials has been established in terms of the traditional 
extremal problem for the weighted logarithmic energy on a compact set of the plane. 

The main goal of this paper is to consider the weak-* asymptotics of the Heine-Stieltjes and 
Van Vleck polynomials in the general setting of ^ C C and pk £ C, which leads to a very different 
electrostatic problem - equilibrium problem in the conducting plane (with a finite exceptional set of 
points). It is essentially known that zeros of Heine-Stieltjes polynomials present a discrete critical 
measure ~ saddle point of the discrete energy functional. A continualization of this notion leads to 
a concept of "continuous" critical measure, i.e. critical point of the usual energy functional defined 
on Borel measures with respect to a certain class of local variations. 

We prove (Section [?]) that weak limit of discrete critical measures is a continuous critical measure 
(as the number of atoms or mass points tends to infinity). Thus, discrete critical measures are limit 
distributions of zeros of the Heine-Stieltjes polynomials. 

To complete the description of the limit zero distributions of these polynomials we have to study 
more deeply the set of continuous critical measures. The problem, rather complex, is connected 
to many other classical problems of analysis, and has potentially a large circle of applications. 
In Section [6] we mention a few connections, in particular, to minimal capacity problem and its 
generalizations . 

In Section [5] we characterize critical measures in terms of trajectories of a (closed) rational 
quadratic differential on the Riemann sphere; for completeness of the reading we summarize basic 
results on quadratic differentials in Section [4] Further investigation of such differentials in carried 
out in Sections [s] (case p = 2) and 9] (general case). 

In the following two sections, [2 and[3| we discuss in some detail the concepts of the discrete and 
continuous equilibrium. 
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2 Discrete and continuous extremal measures 



2.1 Stieltjes electrostatic model: discrete equilibrium 

We denote by ^Mn the class of uniform discrete measures on C, 

^n= \Y,S,„Zk£c\ , and dJl"^ \JdJln, 

lfc=l J n>l 

where 6x is a unit mass (Dirac delta) at x. With any polynomial P{z) = YVj=i{^ ~ Cj) we associate 
its zero counting measure 

n 

i=i 

where the zeros are counted according to their multiplicity. 
For /i = Y12=i ^Ck G 2^ we define its (discrete) energy 



£{fi) = ^log ^ , 



(if two or more Cj's coincide, then = +oo). Additionally, given a real- valued function (external 
field) if, finite at supp(;u), we consider the weighted energy 



n 



£M = ni^)+2Y.^iCk). (2.1) 
fc=i 

In the above mentioned paper [73] Stieltjes introduced the following extremal problem. For 
fixed subset A = {a^, . . . , Op} C M, ag < • • • < fflp, values Pk ^ 0, k = 01, ■ ■ ■ ,p, and an arbitrary 
vector n = (ni, . . . , Up) £ (where = N U {0}), define |n| = ni + • • • + Hp, Aj = [aj^i, aj], 
j = 1, . . . ,p, and A = U^^^^Aj = [ao, Op]. Consider the class of discrete measures 

9Jt|„|(A,n) = {/i G : supp(^) C A, fi{Aj) = nj, j = l,...,p}, (2.2) 
and the external field 

p 

ip{x)=Re{'^{x)), cD(x)=V^log . (2.3) 

^-^ 2 x — aj 

j=0 J 



We seek a measure /i* = n*{n) minimizing the weighted energy ( |2.1[ ) in the class 9Jl|„|(A, n): 

£:^(^*) =min{£:^(^) : ^ G 9Jl|^|(A, n)} . (2.4) 
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In other words, we place rij unit electric charges on the conductor and look for the equilibrium 
position of such a system of charges in the external field ip, if the interaction obeys the logarithmic 
law. 



Stieltjes proved that the global minimum (2.4) provides the only equilibrium position, and that 



the zeros of the solution y = Qn of (1.3) are exactly points of the support of the extremal measure 



H* in (2.4): v{Qn) = Actually, //* provides also the unique component-wise or point-wise 



minimum of £^ ( "Nash- type" equilibrium) . 



The Stieltjes equilibrium problem (2.4) is a constrained one: the constraints are embedded in 
the definition of the class 9Jt|„|(A,n). A classical non-constrained version of the same problem 
leads to the (weighted) Fekete points. Given a compact A C C and n G N, we want to find 
^* G 9Jt„(A) = G 9Jt„ : supp(^) C A} with 

= min{£:^(/i) : G aJt„(A)} . 

Stieltjes' model for the hypergeometric case {p = 1) provides the well known electrostatic in- 
terpretation of the Jacobi polynomials. Zeros of the Jacobi polynomials Pn"'^^ are also weighted 



Fekete points for A = [—1, 1] and <f{x) = log |^^3X| + log |^_)_^| . Similarly, zeros of Laguerre 
and Hermite polynomials are weighted Fekete points for A = [0, -|-oo), (p{x) = log |^ + f and 

2 

A = R, (p{x) = respectively. It was pointed out in [30^ that zeros of general orthogonal polyno- 
mials with respect to a measure on M may be interpreted as weighted Fekete points with an external 
field = in general depending on the degree n. 

Besides its elegance, the electrostatic model just described allows to establish monotonicity 
properties of the zeros of the Heine-Stieltjes polynomials as function of the parameters pk- Fur- 
thermore, the minimization problem for the discrete energy it is based upon, admits substantial 
generalizations (one of them is subject of the present paper). The problem of the limit distribution 
of the discrete extremal points as n — > oo leads to the corresponding continuous energy problems. 

2.2 Extremal problem for Borel measures: continuous equilibrium 

We denote by Ai (resp.. Mm.) the set of all finite positive (resp., real) Borel measures p with compact 
support supp(/i) C C. Hereafter, |^| stands for the total variation of /U G Mr, and \\fi\\ = |^|(C). 
For n G N, let A4n = {fJ. C M : = n} be the set of positive Borel measures with total mass n 
on C. 

With every measure n G A4r we can associate its (continuous) logarithmic energy 

E{p) [ [ log ^^ d^,{x)d^^{y) . (2.5) 
JJ \x-y\ 

Given the external field if G L^(|//|), we consider also the weighted energy 

E^ip,) = E{n) +2 y 99 d/x . (2.6) 



6 



If r is a subset of C, we denote by M.{T) (resp., A4K(r)) the restriction of the corresponding 
famihes to measures supported on T. Again, a standard extremal problem of the potential theory 
is to seek for a global minimizer X-p^ip € Mi{K) such that 



E^{Xr,^) = p = mm{EM- Mi{K)} . (2.7) 

It is well known that under certain conditions on ip this minimizer X-p^ip exists and is unique; it 
is called the equilibrium measure of T in the external field ip, see e.g. [61] for further details. For 
if = 0, measure Ar = Ar,o is also known as the Robin measure of T. 

In terms of the extremal constant p we can also define the weighted (logarithmic) capacity of K, 
cap^(O) = e^f. For ip = we simplify notation writing cap(O) instead of capo(r2). If cap(r2) = 
then O is a polar set. Observe that E(ii) = +oo for any /i G so that any finite set is polar. 

There is a number of properties characterizing the equilibrium measure Ar,!^. For instance, if 
we define the logarithmic potential of /x E Aic by 

U'^iz)'^ I log^-L_dpit), 

then up to a polar subset of K, 

U'^-^{z)+^{z)l=P*^ if.Gsupp(/.^), 
[>/>*, if zGF, 

where p* is a constant related to p and ip. Furthermore, if F and (p are sufficiently regular, 

'u^^-'^{z) + ip{z)) = max min (C/^(z) + 99(2)) . (2.9) 



mm 

2:er V ' ' ' ' V tJ.eM(r), zer 

liMli=l 



This max-min property is a basis for applications of the equilibrium measure in the asymptotic the- 
ory of extremal (in particular, orthogonal) polynomials, see [22], [50], [5^, and also the monograph 

Like for the discrete measures, we will consider general external fields of the form 

(p{z) = Re$(z), 

where $ is analytic, but in general multivalued. What we require in the sequel is that is 
holomorphic in C \ .4, allowing further construction below. 

Further generalizations of this construction can be obtained either considering several measures 
on respective sets interacting according to a certain law (vector equilibrium) |i23j, or including 
additional constraints. For instance, prescribing an upper bound on the density of the extremal 



measure on F in (2.7) we obtain the so-called constrained equilibrium [T7], [5^, relevant for the 



asymptotic description of polynomials of discrete orthogonality. Another way is to impose in (2.7) 
the size of p on each component of F, such as it was done in [45^: if ^ = {gq, . . . , ap} C M, 
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ao < ■ ■ ■ < Up, Tj = [aj_i, Cj], j = 1, . . . ,p, T = U^-^Tj = [uq, Up], and M is the standard simplex 



in RP-i, 



M= le = iei,...,0p) : 0i>O, i = l,...,p, and ^0i = l 



1=1 



then for each = (Oi, . . . , Bp) £ J\f we can consider the global minimum of the weighted energy 
Eip{-) restricted to the class 

Mi{T,e) = {fiGMi: supp(/i) C r, fi{Tj) = 9j, j = 1, . . . ,p - 1} . 



Again, for (f like in (2.3) with pj > there exists a unique minimizing energy, Ar,<^(0). 
2.3 Relation between discrete and continuous equilibria 

The transfinite diameter of a compact set T is defined by the limit process when the number of 
Fekete points tends to infinity. It was Polya who proved the remarkable fact that the transfinite 
diameter of T is equal to its capacity. Fekete observed further that the normalized counting measure 
of Fekete points converges to the equilibrium (Robin) measure of F. For the weighted analogue of 
this result, see [£1, Ch. III]. 

The connection between the discrete and continuous equilibria allowed to use the Stieltjes model 
in [l5] in order to obtain in this situation the limit distribution of zeros of Heine-Stieltjes polyno- 
mials. Namely, if for each vector n = (ni, . . . , Up) G we consider the discrete extremal measure 
p*{n) introduced in (2.4), and assume that |n| — > oo in such a way that each fraction nfc/|n| has a 
limit. 



lim f4 = ^j, i = l,- 



|n|-+oo \n\ 

then /i*(n)/|n| weakly converges to the equilibrium measure Ar,o(^)i T = [oq, a^], 6 = {9i, . . . ,9p), 
defined in the previous section. In a certain sense, this can be regarded as a generalization of the 
just mentioned classical result of Fekete. 



3 Discrete and continuous critical measures 

According to a well-known result of Gauss, there are no stable equilibrium configurations (i.e. local 
minima of the energy) in an a conducting open set under a harmonic external field. Unstable 
equilibria usually do not attract much attention from a point of view of Physics. However, as we 
will show further, they constitute a rich and relevant object that appear naturally in many fields of 
analysis. 

We introduce now the concept that plays the leading role in this paper: the family of measures 
providing saddle points for the logarithmic energy on the plane, with a separate treatment of the 
discrete and continuous cases. 
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3.1 Discrete critical measures 

We start with the fohowing definition: 

Definition 3.1 Let O be a domain on C, A C a subset of zero capacity, and (/? be a real- valued 
function in 17 \ A measure 

n 

A' = E '^Cfc ^ ^ for i + j, (3.1) 

k=l 

is a discrete {A, if) -critical measure in if supp(/x) C C \ .A, and for the weighted discrete energy 
<?<^(/i) = SipiCi, ■■■ Xn) we have 

gradf^(Ci,...,a) = 0, (3.2) 

or equivalently, 

-£:,(Ci,...,.,...Cn)U.=0, k = l,...,n, g-^ = ^[Q^-^Q^ 

More generally, if 99 = Re$, where $ is an analytic (in general, multivalued) function in Q with a 
single- valued derivative then this definition does not need any modification. 

In the sequel we omit the mention to 17 if 17 = C. 



The following proposition is just a reformulation of Eq. (1.3) in this new terminology: 

Proposition 3.2 Assume that A = {ao,ai,. • • ,ap}, p G N, is a set of pairwise distinct points on 
C, and the external field ip is given by ( |2.3[ ), where pk's are the residues in (1.6). Then 

n 

l' = Y.ku^ ^n, C. ^ Ci for % + j, (3.3) 

k=l 

supported on C\A, is a discrete {A, if) -critical measure if and only if there exists a polynomial 
Vn € IPp-1 such that y{z) = yn{z) = YYk=ii^ ~ Cfc) "i-s 0, solution of the differential equation (1.3). 

In other words, discrete {A, (/9)-critical measures with external field generated by complex charges 
fixed at A correspond precisely to zeros of Heine-Stieltjes polynomials. 

Proof. A straightforward computation shows that for z ^ 

2 — log \z — w\ 



dz z — w 

Hence, 

2 » 



,^^f(C.....,C..) = -2Aglog|C.-C.I = -2E^. 
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On the other hand, the multivalued function ip has a single- valued derivative given by (see (|2.3|)) 



E 



Pj 1 



2 Z — Gj 



Thus, using the notation from (1.6), we can rewrite condition (3.2) as 



B 



+ ^(Cfc)=0, k = l,...,n, 



C, A 



(3.4) 



and with y(z) = Yli=i{^ ~ Ci) this identity takes the form 

(^ + f)(a) = 0' k = l,...,n. 

As a consequence, polynomial 

A{z) y"{z) + B{z)y'{z) e P„+p_i , 
is divisible by y, so there exists a polynomial Vn £ IPp-i such that 

A{z)y"{z) + B{z)y'{z) = Vn{z)y{z), 

which concludes the proof. 



(3.5) 



□ 



In the sequel we will make use of the following uniform boundedness of the supports of the 
discrete critical measures, corresponding to a sequence of external fields of the form 



V9n = Re$„, $„(z) = - ^ ^^^y^ log(2; - afc) , 



(3.6) 



fc=0 



where Pk{n) £ C. 

Proposition 3.3 Let /i„ G 971^, n G N, be a discrete {A, (pn)- critical measure corresponding to an 



external field (3.6). // 



liminfRe > ^-^^ > — , 
n ^ n 2 

k=0 



(3.7) 



then |J„ supp(/i„) is bounded in C. 



In other words, if we assume that in (1.3) the coefficient B = Bn may depend on n, but Bn/n is 



bounded (in such a way that (3.7) holds), then the zeros of the Heine- Stieltjes polynomials are also 
uniformly bounded. 
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Proof. Let Hn = Sfc=i ^Cfc(") ^ ^^'^ assume that iCiC?^)! > • • • > |Cn('^)|- Since |Ci(^^)| > 0, by 
(pl), 



But 



so that 



j=2 ^'^^^^ ^ fc=0 



Ci(n) 



IO(n)/Ci(n)|<l 



Ci(n) - afc ■ 



Re 



l-0(n)/Ci(n) 



> 1/2, 



n 



1 ^ 

— Re^pfc(n) 



Ci(n) 



A:=0 



Ci(n) - ttk 



< 



Hence, if Ci('t-) ~^ oo along a subsequence of N, then 

liminfReV^< ' 



k=0 



n 



2' 



which contradicts our assumptions. 



□ 



Remark 3.4 It was proved in [66] that for a fixed of the form ( |3.6| ) (that is, pk{n) = Pk, 
k = 0, . . . ,p), the zeros of the Heine-Stieltjes polynomials accumulate on the convex hull of A. 



Remark 3.5 Condition (3.7) is in general necessary for the assertion of Proposition 3.3 Indeed, 
for p = 0, oo = 0, and 

Tl — 1 

fn{z) = — ^ log|z|, 

any discrete uniform measure supported at the scaled zeros of unity, that is. 



f^n = Yl '^a(n) e ^n, (kin) = Cn e'"''=/", Cn G C \ {0}, 
k=l 



is (^, (/9„)-critical, which is easily established using (3.4) and (3.5). Obviously, for — > oo the 
support of fin is not uniformly bounded in n. 



3.2 Continuous critical measures 

Unlike in the discrete case, we provide now a variational definition for the continuous critical 
measure. 

Any smooth complex-valued function h in the closure $7 of a domain 0, generates a local variation 
ofnhy z ^ = z + th{z), t e C. It is easy to see that z i— > z* is injective for small values of 
the parameter t. The transformation above induces a variation of sets e i— > e* = {z* : z £ e}, and 
(signed) measures: p ^ /Lt*, defined by /U*(e*) = /u(e); in the differential form, the pullback measure 
|U* can be written as (i^*(x*) = dp{x). 
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Definition 3.6 Let O be a domain on C, „4 C O a subset of zero capacity, and 99 be a real-valued 
function in \ We say that a signed measure /U G is a continuous {A, (p)-critical if for 

any h smooth in 17 \ ^ such that h\ ^ =0, 

d 



A 



dt 



t=Q 



(3.8) 



Furthermore, il ip = Re$, where ^> is an analytic (in general, multivalued) function in Q with a 
single- valued derivative then this definition does not need any modification. 

Again, we omit the mention to 17 if 17 = C. In what follows we will always mean by an {A, ip)- 
critical measure the continuous one, satisfying Definition 3.6 Furthermore, in order to simplify 
notation, we speak about an A-critical measure meaning a continuous {A, (/9)-critical measure with 
the external field 93 = 0. Observe that if ^ 7^ 0, this notion is nontrivial. 
A particularly interesting case is treated in the following Lemma: 



Lemma 3.7 // 93 = Re^, and $ is analytic in a simply connected domain O, then condition (3.8) 
is equivalent to 

/^(/i;/i) =0, 

with 



f f Kx) - Hy) 

x-y 



dn{x)dfi{y) 



^'{x) h{x) dfi{x) 



(3.9) 
(3.10) 



Proof. It is sufficient to show that 

E^{fi')-E^{f,)=-Re{tf{fi; h) + 0{t^)} 



We have 
so that 



log 



dfi\x')dfi\y'^ 



\x'- — y 



log 



1 



log 



1+t 



\{x-y)+t{h{x)-h{y))\ 
h{x) - h{y) 



dfi{x)dfi{y) , 



Re 



log 1 + 1 



x-y 

h{x) - h{y) 



x-y 



dfi{x)dn{y) 

dix{x)dii{y) 



On the other hand, 

(/p(x*)d^*(x*) 



ip{x) dfi{x) 



ip{x + th{x)) dfi{x) 



ip{x) d^{x) 
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Taking into account the behavior of log(l + x) for small x, we conclude that as i — > 0, 



E^ifi') - E^i^i) = -Re [t 



h{x) - h{y) 



+ 0{t^) dfi{x)dfi{y) 



x-y 

2 Re I {t<^'{x)h{x) + 0{t^)) dfi{x), 



and the statement follows. 



□ 



Remark 3.8 For a finite set A and the external field given by (2.3), the discrete (/9)-critical 



measures fit into the same variational definition as their continuous counterparts, as long as we 



replace in (3.8) the continuous energy Eip{^) by £ip{^). 



Indeed, arguments similar to those used in the proof of Lemma 3.7 show that for ^ in (3.3), the 
condition 



written for 



yields 



d_ 
di 



0, 



(3.11) 



^(C) 



z ^ A, 



B(z) 

{Q-z){Q-z) ' A{z) 



+ 



fete. 



A{z) ' 



where D is a polynomial. In particular, the residue of the left hand side (as a function of z) is at 
Cfc, k = 1, . . . ,n; setting y{z) = YYi=i(-^ ~ Ci)) arrive again at the system (3.5 ). 



w 



And viceversa, using the chain rule it is easy to show that the condition (3.2) implies (3.11). 



Critical measures constitute an important object; for a finite set A the natural description of 
their structure is in terms of the trajectories of quadratic differentials. In the next section we give 
an abridged introduction to quadratic differentials on the Riemann sphere in the form needed for 
our purposes. For a comprehensive account on this theory see for instance |31l l56| [711 EH]. 



4 Rational quadratic differentials on the Riemann sphere in a nut- 
shell 

Let A and V be monic polynomials of degree p + 1 and p — 1, respectively, with A given b y (|1.2[ ) 
with all Ofc's pairwise distinct. The rational function V/A defines on the Riemann sphere C the 
quadratic differential 

zu{z) = -^^{dzf. (4.1) 
The only singular points of uj (assuming that the zeros of V and A are disjoint) are: 
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• the points S A, where zu has simple poles (critical points of order —1); 

• the zeros of V of order A; > 1, where w has zeros of the same order; 

• the infinity, where w has a double pole (critical point of order —2) with the residue —1. 

The rest of the points in C are the regular points of zu, and their order is 0. All singular points of 
order > — 1 are called finite critical points of zu. 

In a neighborhood of any regular point zq we can introduce a local parameter 




C = e(^)= / I \l~dt, (4.2) 

in terms of which the representation of zu is identically equal to one. This parameter is not uniquely 
determined: any other parameter ^ with this property satisfies ^ = it^ + const. Function ^ is called 
the distinguished or natural parameter near zq. 

Following [56, and [74J, a smooth curve 7 along which 

-V{z)/A{z) {dzf > ^ Im^(z) = const 

is a horizontal arc of the quadratic differential w. More precisely, if 7 given by a parametrization 
z{t), t G (a, /3), then 

-^i<t)) (^) >o, te{a,p). 



A maximal horizontal arc is called a horizontal trajectory (or simply a trajectory) of zu. Analogously, 
trajectories of —w are called orthogonal or vertical trajectories of w; along these curves 

V[z)/A{z) {dz f > 4^ Re^(z) = const . 

Any simply connected domain D not containing singular points of zu and bounded by two 
vertical and two horizontal arcs is called a zu-rectangle. In other words, if ^ is any distinguished 
parameter in D, then £,{D) is a (euclidean) rectangle, and D 1— > £,{D) is a one-to-one conformal 
mapping. Obviously, this definition is consistent with the freedom in the selection of the natural 
parameter ^. 

We can define a conformal invariant metric associated with the quadratic differential zu, given 
by the length element \d^\ = \y^V/A\{z)\dz\; the tu-length of a curve 7 is 

(z) \dz\ ; 

(observe that this definition differs by a normalization constant from the definition 5.3 in [Mj). Fur- 
thermore, if is a simply connected domain not containing singular points of zu, we can introduce 
the ro-distance by 

dist(2;i, Z2;zu,D) = inf{||7||^ : zi,Z2 S 7, 7 C D} . 
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Trajectories and orthogonal trajectories are in fact geodesies (in the tu-metric) connecting any two 
of its points. Indeed, according to [56, Thm. 8.4], in any simply connected domain D not containing 
singular points of w, a trajectory arc 7 joining zi with Z2 is the shortest: if Li, L2 are the orthogonal 
trajectories through zi and Z2, respectively, then any rectifiable curve 7 that connects Li with L2 
in D satisfies 

II7IU < II7IU • 



Figure 1: The local trajectory structure near a simple zero (left) or a simple pole. 

The local structure of the trajectories is well known (see the references cited at the end of the 
previous Section). For instance, at any regular point trajectories look locally as simple analytic 
arcs passing through this point, and through every regular point of lu passes a uniquely determined 
horizontal and uniquely determined vertical trajectory of w, that are locally orthogonal at this 
point \74:\ Theorem 5.5]. If z is a finite critical point of tu of order k > —1, then from z emanate 
k + 2 trajectories under equal angles 2TT/{k + 2) (see Figure [T|. 

In the case of a double pole, the trajectories have either the radial, the circular or the spiral 
form, depending whether the residue at this point is negative, positive or non-real, see Figure [2} 
In particular, with the assumptions on A and V above all trajectories of the quadratic differential 



(4.1) in a neighborhood of infinity are topologically identical to circles. 

The global structure of the trajectories is much less clear. The trajectories and orthogonal 
trajectories of a given differential w produce a transversal foliation of the Riemann sphere C. The 
main source of troubles is the existence of the so-called recurrent trajectories, whose closure may 
have a non-zero plane Lebesgue measure. We refer the reader to [71] for further details. 

A trajectory 7 is critical or short if it joints two (not necessarily different) finite critical points 
of w. The set of critical trajectories of w together with their endpoints (critical points of w) is the 
critical graph of w. Critical and closed trajectories are the only trajectories of w with finite w- 
length. The quadratic differential w is called closed if all its trajectories are either critical or closed 



15 



Figure 2: The local trajectory structure near a double pole with a negative (left), positive (center) 
or non-real residue. 



(i.e. all its trajectories have a finite tu-length). In this case the trajectories of w that constitute 
closed Jordan curves cover the whole plane, except a set of critical trajectories of a plane Lebesgue 
measure zero; see e.g. Figure [3] for a typical structure of such trajectories. 



If the quadratic differential ( |4.1[ ) with A given by ( 1.2 ) is closed, there exists a set F of at most p 
critical trajectories of w such that the complement to F is connected, and \fVjA has a single-valued 
branch in C \ F. 

5 Critical measures in the field of a finite system of fixed charges 

In what follows we fix the set of p -|- 1 distinct points A = {ao; • • • i Op} C C and consider the basic 
domain O = C \ ^ = {ag, oi, . . . , Op}, and an external field Lp of the form 

fc=0 k=\ ^ ^ ' 



where we have used notation from (1.6). If {poi • • • i Pp\ C M, then this external field corresponds to 



the potential of a discrete signed measure supported on A: 



^{z) = U'^{z), a = ^^5,, G2)Tp+i. (5.2) 

k=0 



However, if any pk G C \ M, then ip is not single-valued in C \ nevertheless, the notion of an 



{A, 93)-critical measures for this case has been discussed in Definition 3.6 In particular. Lemma 3.7 
applies. 

In this section we state and prove the main structural theorem for (A, (/3)-critical measures, which 
asserts that the support of any such a measure is a union of analytic curves made of trajectories of a 
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rational quadratic differential. On each arc of its support the measure has an analytic density with 
respect to the arc-length measure. Finally, we describe the Cauchy transform and the logarithmic 
potential of an {A, (/?)-critical measure. 



5.1 The main theorem 



According to (5.1), A is exactly the set of singularities of the external field (p, except for the case 
when pk = for some k £ {0, . . . ,p}. In such a case we do not drop the corresponding from the 
set A; it remains as a fixed point of the class of variations (Definition 3.6). However, the status of 
the point £ A with = is different from the case pk ^ 0, see next theorem. 



Theorem 5.1 Let A = {oq, ai, • . . , flp} and (/? given by (5.1). Then for any continuous {A,ip)- 
critical measure p, there exists a rational function R with poles at A and normalized by 

i?(z) = (-)Vo(4)' ^ = MC) + ^X^P,-, (5.3) 

such that the support supp(p) consists of a finite union of trajectories of the quadratic differential 
w{z) = —R{z)dz'^ . If all pj £ M, then vj is closed. 

If in the representation (5.1), pj = 0, j £ {0,1,...,^}, then aj is either a simple pole or a 
regular point of R; otherwise R has a double pole at aj. 

The proof of this theorem reduces to two lemmas below. The first of them deals with the 
principal value of the Cauchy transform 



C^(z) = lim / ^— dp{x) 

of the {A, (/3)-critical measure p. 

Lemma 5.2 For any {A, <p)- critical measure p there exists a rational function R with properties 
listed in Theorem 15. Jl such that 

{C>'{z) + '^'{z)f = R{z) mes2-a.e., (5.4) 

where mes2 is the plane Lebesgue measure on C. 



1 



Proof. Assume that p is an {A, (/3)-critical measure for ip like in (5.1 ). We will actually show that 
(5.4) is valid at any point z £ C where the integral defining is absolutely convergent. It is well 
known that at such a z, 



lim / M^ = o, (5.5) 

r^0+ 7|^_^|<,, \X - Z\ 

and this property holds a.e. with respect to mes2- 
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At this point we would like to emphasize that the variational arguments we present next would 
be significantly simpler for points z € supp(/i); in this case one can use the standard Schiffer 
variations = — z) (see e.g. the Appendix of [H]), as it was done in the original paper 

|59j (for the logarithmic potentials), and subsequently in |33| Chapter 8] and pi] (for the Green 
potentials). In the present situation we do not have any a priori information about supp(/u). In 
order to address the problem of a possible pole on the support of ^u, we modify the variation by 
making h{C^ = in a neighborhood of z (see (5.7)). The function 6{C) introduced below is meant 
to preserve the smoothness of h. 

Fix z G C satisfying (5.5), and for r > denote Dr = {C G C : |C — z| < r}. Function 



m(r 1 



^ (2?r) is continuous from the left and monotonically increasing, so that the subset 



A 



II \ ^ m(r + e) — m(r — e) 
r G (0, 1) : m (r) = lim exists 



has the linear Lebesgue measure 1. 

For r G A and e G (0, 1) define the "smooth step" function 



A(a:,e) 



0, 

(x-l-e)2(x-l + 2e) 



4e3 



.1, 



if < X < 1 - e, 

, ifl — e<x<l+e, 
if X > 1 + e. 



It is easy to see that A(-,e) G C^(]R+) and that |^ A(x,e)| < 1/e for all e > 0. Using this function 
we define on C the function 



^"Af'^,£,, 



and consider the condition (3.9) with the following particular choice of h: 

h{C)=he{C,r) = ^e{C,r,e). 



For the sake of brevity we use the notation 



SO that T>r{i-e): ^r,e and J-r,e provide a partition of C. Furthermore, by construction 



(5.6) 



0, if CG2?r{i-.), 

^(0 = <! A{0 



(5.7) 
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Consider first 

he{x;r) - he{y;r) 
x-y 



dll{x)d^{y) = X Pr(l-e)) + I{^r,e X ^r,e) + I{^r,e X ^r,e) 



2/(P,(i_,) X ICr,e) + 2/(P,(i_,) X J^r,e) + 2/(/C,,, X ^,,,), 

where means the integral in the l.h.s. taken over the set Observe that by ( |5.7[ ), /(Pr(i-e) x 

^?r(l-.)) = 0. 

Let C £ l^r,e] since 
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6'(C) = - A' 



1 ./^IC-^ 



d 



r J dC r 



1 ^>f\C-^\ A C-^ 



IC-^I 



we have 



|grade(C)|| 



d 



In consequence, for x,y E JCre, 



he{x;r) - he{y;r) 



0{C) 



< 



re 



x-y 



< 



const 



re 



(5.8) 



where the constant in the right hand side is independent of e. Obviously, by definition of h we have 
that this inequality is valid (with a different constant) if x G ICr^s and y lies on a compact subset of 
C. 



From (5.8) we conclude that 



\I{)Cr,e X )Cr,e)\ 



K{x;r) - he{y;r) 



x-y 



dii{x)dix{y) 



Taking into account that r E A, we have that 

hm 

s^o+ e 



^ const 2 



(5.9) 



2r m'{r), 



(5.10) 



so by ([5^, /(/Cr,£ X ICr,s) = o(l) as e ^ 0+. 
Consider now x G fCr^e a-^d y G P^.e- Then 

he{x]r) - he{y;r) ^ he{x;r) ^ 

x-y x-y (x — z){x — y) 

Consider two cases. If \y — z\ < r(l — 2e), then 



A{x) 



9{x) 



he{x;r) - he{y;r) 



x-y 



< 



const 



< 



const 



r{l-e)\x-y\ r{l - e){\x\ - \y\)' 
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Hence, with a different constant, 



he{x;r) - he{y;r) 

cG/C,.,e, |y-2|<r{l-2e) X-y 



dfi{x)dn{y) 



< 



const 



r{l+£) /■r(l-2e) 



r(l~e) Jo 



St 
t - S 



dsdt: 



the double integral in the r.h.s. is explicit, and it is straightforward to verify that it is o(l) a e ^ 0+. 

If on the contrary r(l — 2e) < |y — 2;| < r(l — e), with x S /Cr,e we can use the estimate (5.10), 
which yields 



hs{x;r) - he{y;r) 

€Kr,e, r{l-2e)<\y-z\<r{l~e) X — y 



dfi{x)dfi{y) 



re 



and again by (5.10), the r.h.s is o(l) a e ^ 0+. Gathering the last two estimates we conclude that 

X }Cr,e) = o(l), as e ^ + . 



Similar considerations are obviously valid for I{fCr^e ^ ^r.e)- 
Summarizing, 

^hm^ (/(/C,,e X lCr,e) + 2/(2?r(l-e) X /C,,,) + 2I(/C,,, X :Fr,e)) = 0, 



(5.11) 



and we conclude that 

he{x;r) - he{y]r) 



lim 



x-y 



dfi{x)dfi{y) 



+ 2 



ho{x;r) - ho{y;r) 

\x—z\>r, 2|>r X y 

ho{x;r) - ho{y;r) 



where 



' \x— z\>r,\y—z\<r X y 

/i(r)+2/2(r), 

0, if IC - -^1 < r, 

A{0 



dfi{x)dfi{y) 



dn{x)dij,{y) 



(5.12) 



-, if IC — -^1 > T. 



ho{C;r) = < 

Let us analyze the behavior of each integral as r — > O"*" separately. First, 

Aix) 



h{r) 



\x-z\>r,\y-z\<r {x z){x y) 



dfi{x)dn{y). 



Observe that by (5.5), 



\x~z\>r 



A{x) 



X — z 



\y—z\<r \x ~ y\ 



dfi{x) < +00, 



20 



so that applying Fubini's theorem we conclude that 

A{x) ( f dniy) 



.ir)^f ^ f 

J\x-z\>r X — Z \J\, 



Using again (5.5) we obtain that 



Ix— 2|>r Z \J \y~z\<r V 



lim l2{f) = 0. 



dfi{x) . 



On the other hand, by (|5.7|) 
/i(r) = 



A{x) 



My) 



dfi{x)d^i{y). 



\x-z\>r,\y-z\>r \{x - z){x - y) {y - z){x - y) _ 
The identity 

A{x){y - z) - A{y){x - z) + A{z){x - y) = {x - y){x - z){y - z)D{x, y, z) 
is immediate, where 

D{x, y, z) = aQ{x, y) + ai{x, y)z -\ h ap-2(a:, y)zP~'^ + z^"^ 

is a polynomial of degree < p — 1 in each variable. Hence, 

A{x) A{y) Aiz) 



{x-z){x-y) {y-z){x-y) {x - z){y - z) 



+ D{x,y,z) . 



Using it in (5.13) we get that 



where 



Thus, by (5.12), 



lim lim 

■»0+,r6Ae^0+ 



hm h{r) = D,{z) - A{z) (C^(z))^ 



Di{z) = j j D{x,y,z) dn{x)dn{y) 
he{x]r) - he{y;r) 



x-y 

In a similar fashion we can analyze 



d^Ji{x)d^Ji{y) = Di{z) - A{z) {C^{z)y . 



^'{x)he{x;r) d^i{x) = [ / + / ] ^'{x)he{x;r) diJ.{x). 
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Again estimates on lCr,£ and (5.10) show that 



lim / ^'{x)he{x;r) dfj,{x) = / ^>'( 

-^O+J J\x-z\>j- 



A{x) 



\x—z\2 

Taking into account (|5.1l) we can rewrite the right hand side as 



x — z 



dfj,{x). 



B(x) 



x — z 



dfi{x) 



1 



B{x)-B{z) ^ 1 

^ J\x-z\>r X — Z Z 



\x—z\>r 



X — Z 



dfi{x). 



Thus, 



where 



lim lim / ^' {x)hs{x;r) dii{x) 



{D2{z) + B{z)C^^{z)), 



D2{z) 



B{x) - B{z) 



X — z 



dfi{x) 



(5.18) 



is a polynomial of degree < p — 1; it is = 0if99 = 0. 



Combining this last identity with (5.17) and using (3.10) we get that 

lim lim f^ifi- h) = Di{z) - A{z) + D2{z) + B{z)C^{z). 

r— >0 £— >0+ 



(5.19) 



Since (3.9) is valid for each e > and r € A, we obtain that the right hand side in (5.19) is 0. 
In consequence, 

{A{z)C^'{z)f - A{z)B{z) C^(z) + B\z)/A = A{z) {Di{z) + D2{z)) + B''/A{z) . 
Taking into account ( |5.1[ ) we rewrite this condition as ( |5.4[ ), with 

dcf Di{z) + D2{z) 



R{z) 



Aiz) 



and clearly R has a double pole at z = aj if and only if pj ^ 0. Finally, the normalization condition 
(5.3) follows from considering (5.4) as z — > cxd. This establishes the assertion of the theorem. □ 

The next proposition is a slightly modified version of [51 Lemma 4] by T. Bergkvist and 
H. RuUgard. The original lemma considers only positive measures /j, for which (C^)*^, for certain 
A: G N, is a reciprocal of a polynomial. 

Lemma 5.3 Assume that p, G Mm. is a finite signed Borel measure on the plane whose Cauchy 
transform is such that that there exist rational functions r and R with possible poles at A 
satisfying 



{C + rYiz) = R{z) 



mes2 — a.e. 



(5.20) 
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Then fj, is supported on a union of analytic arcs, that are trajectories of the quadratic differential 
w{z) = —R{z)dz^ , with possible mass points at A. 

If all poles of r are simple and have real residues, then additionally w is a closed differential, 
and the number of connected components of supp(/i) is finite. 

Finally, if is a positive Borel measure, then any w-rectangle not containing the zeros or poles 
of R intersects supp(^) along at most one analytic arc. 

Proof. For the quadratic differential w consider a tu-rectangle D (see the definition in Section |4]), 
disjoint with A and not containing the zeros of R. We select in this rectangle a holomorphic branch 
of and a distinguished parameter 

i = i{z) = ^/mt)dt, (5.21) 

which is a conformal mapping of L> onto D = £,{D). 
Let us define in D the following function: 

X(0=sgn(^^(z(e))). 



Hence, x takes only two values, ±1, and for z £ D, 

{C^^ + r){z) = x{az))VR{z). (5.22) 
We have that in the sense of generalized derivatives, for z £ D, 

d d I f d d 



Differentiating in (5.22) and using the chain rule, we get 



Mz) = I (x(e(^))Vi?(.)) = ^ vr{z) = ^ Vii(.). 

Taking into account the definition of x in (5.21) we conclude that z £ D and ^ = ^{z), then 

In particular, the (generalized) partial derivative of x{€) along the vertical axis is zero; \i ^ = x + iy, 
this implies that 

x{Ci = 9{x) X i(y), 

where l(y) is the identity generalized function in the single variable y. Thus, the set 
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is a union of vertical arcs in the .^-plane. It means that the image of the support of /i in D by (5.21 ) 



is made of vertical lines, that is, supp(/i) n D is a union of horizontal trajectories of -cu. Moreover, 



if fi is positive, we get by (5.23) that 



dm 



> in £>, 



so that xiO changes sign at most once in D. In other words, supp(/x) n D contains at most one 
single analytic arc, which is a horizontal trajectory of zu. 
Finally, if all residues of r are real, then 

/z i-z 
+ r) {t)dt = Uf'iz) + Re / r{t)dt + const 

is harmonic and single-valued in C \ .4, which means that the trajectories of tjj are either critical or 
level curves of a harmonic function. Thus, zu is closed. □ 

Remark 5.4 See Figure [3] for an illustration of the statement about the number of intersections 
of a ro- rectangle with supp(/i). Obviously, in this assertion we can replace the ro-rectangle by any 
simply-connected domain that is mapped one-to-one by ^ onto a convex set. 





Figure 3: ro-rect angles intersecting the support of a positive (left) and sign-changing measure 



(right); for further details, see Section 8.4 
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Remark 5.5 Observe that in Lemma 5.3 we do not assume a priori that Eip{ii) < oo, so that mass 
points of fx are allowed. 



An immediate consequence of Lemma 5.2 and Lemma 5.3 



IS 



Corollary 5.6 // in (5.1 ), {po, . . . , pp} C M, then in any connected component of C\ supp(/i) we 
can select a single-valued branch of y/R such that there the formula 



Wiz) + ip{z) 



Re J ^R{t) dt 



(5.24) 



holds. 



5.2 Critical and reflectionless measures 



By Theorem 5.1 supp(/i) is a union of analytic arcs. 



Definition 5.7 We call a point z S supp(/i) regular if there exists a simply connected open neigh- 
borhood B oi z such that B n supp(;u) is a Jordan arc. 

Lemma 5.8 Let fi be an {A, (p)-critical measure. Then the principal value of the Cauchy transform 
of p satisfies 

Ciz) + '^'{z) =0, ze supp(^) \ A. (5.25) 

On any simple subarc of supp(//) measure p is absolutely continuous with respect to the arc-length 
measure, and its density is given by 



dp{z) 



1 



vr 



\/R{z) dz 



(5.26) 



Proof. Assume first that z G supp(^) is regular, and let i? be a simply connected open neighbor- 



hood B oi z such that B n supp(/i) is an open analytic arc not containing A. By Lemma 5.2 



{C^{C) + ^'{Of = R{C), C e C \ supp(^), 

where the Cauchy transform is understood in the strong sense (ordinary integral), so that 

cic) = -^'(c) + Vm (5.27) 

in each connected component of i? \ supp(/x), with an appropriate selection of the branch of the 
square root. 

At a regular point z the boundary values of from both sides of supp(/i). 



C^{z) hm C'^(C), 

C— >2i, CSC\supp(m) 
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are well defined and satisfy the Sokhotsky-Plemelj relations, 



(5.28) 



By (5.27), if the branch of V R coincides on both sides of i? Pi supp(/x), then jj 
impossible. Hence, with an appropriate selection of the branch of in B, 



there, which is 



C^'iz) = -<^'{z)±^/R(zj, z G -Bnsupp(/x), 



and both (5.25) and (5.26) follow from (5.28). 



Finally, if z £ supp(//) \^ is not regular, it must coincide with a zero of the rational function R 



at z at least as the square root of (C 



in (|5.4|). Taking into account the expression for the density of the measure we see that it vanishes 

z). Then, at this point (5.5) holds, as well as formula (5.4). 

□ 



This concludes the proof of (5.25). 



Formula ( 5.25 ) is a direct continuous analogue of property ( 3.4 ) for the discrete critical measures, 

outside 



and it may be proved directly (independently from Theorem 5.1 ) using local variations {h 
of a small neighborhood of the singularity at z 



The original proof of Theorem |5.1| by the second 
author (unpublished) was based on a combination of Lemmas 5.2 and 5.8 The technique from [5j 
used in Lemma |5.3| above streamlines the arguments. 

Observe that for 93 = we obtain from Lemma 5.8 that for any ^-critical measure /i, and for 
any regular point z £ supp(^), C'^{z) = 0, so that = fi-a.e. Measures with this property are 
called reflectionless; see |l9], where they are treated in the context of the geometric function theory. 



Remark 5.9 For measures supported on the unit circle T = {z G C : |z| = 1} we can find in 
literature an alternative definition of the "reflectionlessness" (see [19], [20], [21]), characterized by 
vanishing of the sum of the boundary values of the Caratheodory function of fj, and not of its Cauchy 
transform. For instance, the Lebesgue measure mesi on T is reflectionless in the sense of |l9j, but 
does not satisfy v. p. C"^'^'^^ = on T. In this paper we give the "reflectionless measure" the meaning 
specified above. 

Reflectionless measures on M have their origin in the spectral theory; they are in fact spectral 
measures of Schrodinger self-adjoint operators with reflectionless potentials (see e.g. [12]) and of 
reflectionless Jacobi operators |77] . 

It is immediate to show that there are reflectionless measures with infinite energy. So, the 
following conjecture seems natural: 



Conjecture 5.10 For any positive reflectionless measure with a finite energy we can find a polar 
set A on C such that /x is A- critical. 

A weaker statement is the following. Assume that /i is a reflectionless measure supported on a finite 
set of analytic arcs F = FiU---UFfc ofC, and that ^ is absolutely continuous with respect to the 
arc-length measure. Is it true that this measure is „4-critical, for a suitable finite set A? 
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5.3 Equilibrium conditions in terms of potentials, and the ^-property 

The variational requirements defining a continuous {A, 93)-critical measure /i impose equilibrium 
conditions that we discuss next. We will see that the gradient of the total potential (that is, force) 
vanishes at any regular point of the support of the measure located in the conducting part of the 
plane. However, grad(C/^ + ip) is not continuous across any arc in supp(/i), and we have to consider 
separately the force acting on an element of charge from either side of supp(/i). This leads to 



equality of the normal derivatives, the so called S-property, see (5.29) below. 

Lemma 5.11 The total potential of an {A, f) -critical measure /i satisfies the following properties: 

(i) i/supp(/i) = Ti U • • • U Tk, where Tj are the connected components o/supp(/x), then 

U^{z) + ip{z) = Wj = const, z E Tj, j = 1, . . . , k. 

(ii) at any regular point z £ supp(/i), 

^ (U^ + ^){z) = ^(U^ + ^){z), (5.29) 



5n+ dn 

where n± are the normal vectors to supp(/i) at z pointing in the opposite directions. 
Additionally, if z £ supp(^) \A is not regular, then 

grad(C/^(z) + (/7(z)) = 0. (5.30) 

Reciprocally, assume that a finite real measure fi, whose support supp(/u) consists of a union of 
a finite set of analytic arcs, supp(^) = Ti U • • • U Tyt, satisfies conditions (i) and (ii) above. Then fi 
is {A, (p)-critical. 



Proof. Let /i be an (^, (/?)-critical measure. From Theorem 5.1 and Lemma 5.8 it follows that ^ 



has an analytic density on the regular points of its support, made of analytic curves. Hence, [/^ is 



continuous up to the boundary, and (i) is a direct consequence of Corollary 5.6 and the fact that /u 
lives on trajectories of the quadratic differential —R{z){dz)'^ . 
For any z G C \ supp(/x) we have that 

d 1 



^_ (U^{z) + ^{z)) = -{C>^{z) + ^'{z)), 
and this relation is inherited by the limit values on the Caratheodory boundary of C \ supp(/x). 



Using (5.25) we conclude that on regular points of supp(//) \ A, 



^ {U^\z) + v.(z))+ + ^ {U^{z) + ^{z))_ = 0. 
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Observe that for a real valued function u, du{z) / (dz) coincides, up to a factor 1/2, with the 
(complex) gradient of u. From (i) it follows that grad(C/'^ + ip) is normal to supp(/i), so that 



?rad([/'^(z) + (^(z) 



d 



dn± 



{U^{z) + ^{z)) 



z £ supp(//) \ A, 



and (5.29) follows from the last two equalities. 



Finally, the only possible non-regular points of supp(/i) \^ are the zeros of R, where the density 



vanishes, and we conclude (5.30). 



Let us prove the reciprocal. Assume that ^ G A^k has a bounded support comprised of a finite 
number of smooth linear connected components Tj. Fix a regular point C € supp(//) (without loss 
of generality, z £ Tj), and let again i? be a simply connected open neighborhood of ( such that 
B n supp(;u) is a Jordan arc. It splits B into two disjoint domains, that we denote by B^, so that 



B \ supp(/i) = B+U B-. It follows from (i) and (5.29) that 



U{z) 



'U^'{z) + ip{z)-Wj, 

-{U^^{z) + ip{z)-w,) 
10, 



if z G B+, 
if z G B-, 

ii z £ B n supp(^)-|- 



is harmonic in B. Equivalently, C{z) = dU{z)/{dz) is holomorphic in B. But 

C{z) 



C^(z) + $'(z), 
- (C^(z) + $'(z)) 



if z G B+, 
if z G B-, 



(5.31) 



is continuous in B. It implies that R{z) = {Cfi + <I>')^(z) is holomorphic in B, and in consequence, 
R is analytic at any regular point of supp(/i). Since R is obviously analytic also in C \ supp(/i), we 
conclude that it is in fact holomorphic in C, except for the set of irregular points of supp(;u), that is. 



the endpoints of the arcs comprising supp(yLt). By (5.30) and (5.31 ), R vanishes also at the irregular 



points of supp(^) \ A and at infinity. Since n is a finite measure, has a sub-polar growth at any 
point of C. Altogether it means that i? is a rational function with possible poles at A (of order 
< 2). By Lemma 5.3, measure fi is (^, (/9)-critical. □ 



5.4 Correspondence of critical measures with closed quadratic differentials 



We begin with some general remarks on {A, 93)-critical measures for the case (5.2), when all pfc G M 
and the external field corresponds to the potential of a discrete signed real measure supported on 



A. According to Theorem 5.1, for any (^, (/?)-critical measure fi there exists a closed differential 
-R{z)dz^ such that supp(/x) consists of a finite union of its trajectories. 



This is not a one-to-one correspondence: even in the class ^A the same quadratic differential 
may correspond to a whole family of critical measures. 
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Example 5.12 Let p = 0, oq = 0, and (p{z) = ^ log \ z\ (generated by a charge —1/2 at the origin). 
Then for any r > the normahzed angular (Lebesgue) measure mesi living on the circle \z\ = r is 
{A, (^)-critical. Each such a measure is supported on a trajectory of the same quadratic differential 



(dz) / z (for a discrete analogue of this statement, see Remark 3.5). 



Example 5.13 In a more general situation we can consider an ^-critical measure for an arbi- 
trary configuration A without external field (ip = 0); the trajectories of the associated quadratic 
differential w near infinity are closed Jordan curves. Select any such a trajectory (3 containing in 
the bounded component of its complement both A and supp(^), and denote by /2 the balayage of fi 
onto 13 (see the definition in j6H §11.4]). Then fii = 2'jl — n is another ^-critical measure with the 
same total mass than and such that (3 C supp(/xi); observe that also corresponds to the same 
quadratic differential. The verification of this assertion is a simple exercise. 

What these very basic examples have in common is that in each case when we were able to construct 
more than one ^-critical measure associated with the same quadratic differential, closed trajectories 
in the support of measures were present. Moreover, we have seen that an infinite family of critical 
measures may correspond to the same quadratic differential. This is not the case if we restrict 
ourselves to critical measures with a connected complement to the support. However, even in this 
situation the quadratic differential can give rise to more than one (signed) critical measure, as the 
following example shows. 

Example 5.14 Consider the quadratic diffferential w with 4 simple poles a^^ . . . at the vertices 
of a rectangle, and two simple zeros fi,f2, situated symmetrically at the midpoints of the longest 
sides of this rectangle. The trajectories of such a differential are depicted in Figure |4j We can 
associate to w three different ^-critical measures with a connected complement of their supports. 
Indeed, although the critical trajectories joining poles will always belong to the support of any such 
a critical measure, for the third component of the support we can choose any of critical trajectories 
connecting both zeros Vj. 

All these examples illustrate the general difficulty of the analysis of the correspondence between 
quadratic differentials and critical measures. Nevertheless, we will show below that in the class 
of positive {A, (/?)-critical measures ^ corresponding to ip generated by a positive measure, the 



mapping associating to such a /x the quadratic differential w described in Theorem 5.1 is an injection. 
Moreover, C \ supp(/x) is connected. The assumption of positivity of the mass giving rise to (/? is 
necessary, as the following example (first considered by Teichmiiller in his "Habilitationsschrift" 
r6l) shows. 



Example 5.15 Let A = {0,1}, a = — a(5o, < a < 1 (negative "attracting" charge). Then 
there exists a unique positive critical measure ^ with /u(C) = 1 + a. The corresponding quadratic 
differential is 

_{z-c){dzf^ c = c(a)G(0,l). (5.32) 
z'^[z — 1) 
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Figure 4: Trajectories of the quadratic differential w described in Example 5.14 Bold lines represent 
the critical graph of w. 



The supp(//) is the whole critical graph of this differential, which consists of the segment [c, 1] and a 
closing loop passing through c and enclosing the origin (see Figure [5]). Thus, we do not have critical 
measures with a connected complement to supp(/u). Moreover, if we now consider an external field 
of an opposite sign, ip = U^"^ (positive "repelling" charge), then the corresponding (A,(p)- critical 



measure Jl is associated with the same quadratic differential (5.32), but now supp(/i) = [c, 1], and 
pt(C) = 1 — a. We point out that this is not a mere artificial example: the whole variety of the 
critical measures in a similar situation appears in the asymptotic analysis of the Jacobi polynomials 
with varying non-standard parameters, see e.g. [35] and 



We present here a lemma that will allow us to isolate the cases that are of our interest. 

Lemma 5.16 Let n G TWk he an {A, (p) -critical measure for = U'^ , a £ Mr, and assume that j3 
is a closed contour contained entirely in supp(/i), delimiting the bounded domain $7. Then 

-/i(/3) = 2(/i + a)(17). 

Proof. Let n_ be the unit outer normal vector to (3, and n_|_ = — n_. By Gauss theorem (see 
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Figure 5: Critical graph of the quadratic differential given in (5.32). 

e.g. Theorem 1.1, §11.1 of [61j), 

(f/^ + ^) {z)\dz\ = ^ / U^^+%z)\dz\ = + U 17), 
{U^ + cp) {z)\dz\ = ^ jf ^ U'^+^{z)\dz\ = -(/. + am. 
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By the S'-property (5.29), both integrals in the left hand side are equal, and the Lemma follows 



since (t(/3) = 0. □ 

Proposition 5.17 If <p = , with u > 0, then 

(i) the support of any positive (A, if) -critical measure has a connected complement; 

(a) the correspondence between the positive {A, f)- critical measures and the associated quadratic 
differentials is injective. 



Proof. Statement (i) is an obvious consequence of Lemma 5.16 



Let now be a positive (^, (/?)-critical measure and let tu = —R{z)dz be the quadratic differ- 



ential whose critical trajectories support /i (see Theorem 5.1). We have to prove that w determines 
fi uniquely. 

Let r be the critical graph of w, and denote by Crit{w) the class of all (signed) {A, (/3)-critical 
measures /x corresponding to w and such that C \ supp(/I) is connected; by (i)., fj, S Crit{w). If 
both sides of a critical trajectory 7 C T belong to the boundary of the same connected component 
of C \ r, then either one of these disjoint possibilities holds: 

• 7 belongs to the support of every measure from Crit{w)] 
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7 is not contained in the support of any measure from Crit{w) 



Indeed, by (5.20), any measure Jjl € Crit{w) is recovered from its support using the Sokhotsky- 
Plemelj formulas. Since both sides of 7 belong to a connected complement of C\r, either possibility 
above is determined by the analytic continuation of ^/R: 7 G supp(/i) if and only if ^fR will have 
opposite signs on both sides of 7. This is obviously the case if any closed curve, contained in this 
connected component and joining both sides of 7, encloses an odd number of singular points of w. 

As a corollary, we conclude that any critical trajectory emanating from a simple pole of ro, must 
belong to the support of any Jl £ Crit{w). 

Let now 7 be a critical trajectory joining two different zeros of w. Each side of 7 is a boundary 
of a ring domain filled with closed trajectories of vj. If both sides of 7 are in the boundary of the 
same ring domain, then considerations above apply. So, it remains to consider the case when 7 is 
in the outer boundary of a ring domain J7 (see e.g. the middle arc joining vi and V2 in Figure |4]). 
Let 7 be a closed trajectory from J7, and denote by A the restriction of ;U + cr contained inside 7. 
Since A is by assumption a positive measure, the gradient (flux) of the potential on 7 is directed 
inwards. By continuity, this also happens on 7, so the restriction of /x + cr to 7 cannot be positive. 
In conclusion, 7 does not belong to supp(/u), and it finishes the proof of (ii). □ 

Remark 5.18 In fact, the following stronger uniqueness property holds (formulated in the notation 
introduced in the proof above): if for a quadratic differential TjJ, the set Crit{w) contains more than 
one measure, then none of the measures in Crit{w) is positive. This is the case, for instance, of 



Example 5.14 



A combination of Lemma 5.11 and Proposition 5.17| yields the following addendum to Lemma 5.3 



which has an independent interest: 



Proposition 5.19 Assume that fi £ M is a finite positive Borel measure on the plane whose 
Cauchy transform is such that that 

{C^f{z) = R{z) mes2-a.e. 

for a rational function R. Denote by A the set of poles of R. Then /i is an A-critical measure, and 
it is uniquely determined by R. 



6 Critical measures and extremal problems 

Critical measures are connected in an essential way with a class of extremal problems that lies on a 
crossroad of the geometric function theory, approximation theory, potential theory and some other 
topics that will be mentioned below. We start with one of the oldest problems of that kind. 
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6.1 Chebotarev's continuum 



For a finite set A = {ao, . . . , Up} in C we are interested in the continuum of minimal capacity 
containing A. More precisely, if we denote by J- the family of all continua F C C with ^ C -F, we 
seek T £ such that 

cap(r) = mincap(F). (6-1) 

This problem was raised by Chebotarev (alternative spelling, Tchebotarov) in a letter to Polya 
(see |55]). Grotzsch [26 and Lavrentiev (or Lavrentieff) |39| HQ] proved that there exists a unique 
r* = r*(^) satisfying ( |6.1[ ), and that T* is a union of critical trajectories of a rational quadratic 
differential -R{z){dzf, where R = V*/A, with A given by ((Of, and V*{z) = U'jZliz - v*). We 
call this r* the Chebotarev's compact or Chebotarev's continuum. 

The most recent account on the background of this problem and some of its applications can be 
found in |53| . Here we describe briefly (and inductively) some basic facts about the geometry of T* 
that will be useful in the sequel. 

The case p = 1 (that is, when A contains 2 points) is trivial: in this situation T* is the segment 
[ao,ai] joining both points. 

For p = 2 (the first non-trivial case), V*{z) = z — v* . If ao, ai and 02 are collinear, say 
0-2 £ [oO) then v* = 02, and we are in the situation of p = 1 considered above. Otherwise r*(^) 
is made of three analytic arcs, each emanating from the respective pole aj (j = 0, 1,2) of R, and 
all merging at v* (see Figure |6]). Point v* = v*{A), as function of A, is uniquely defined, but its 
analytic representation is not known (not to speak of the zeros of V* for p > 3), see |38]. In the 
totally symmetric case when Uj lie at the vertices of the equilateral triangle, point v* is just its 
center, and the three curves coincide with the bisectors joining it with the vertices. For further 
geometric properties of F* and v* see |38| Ch. 1]. 

For p = 3 (four poles), in a generic situation the zeros of V* are simple, and the Chebotarev's 
compact consists of 5 arcs; if these zeros coalesce forming a double zero, then r*(^) is made of 4 
arcs (see Figure . When one or both zeros Vj coincide with a pole from A we are left in one of 
the cases previoiisly considered. 

The number of "degenerate" situations growth fast with p, so in the sequel we restrict our 
attention to a generic one, when the points from A are in a general position. This notion of 
genericity (as opposed to some more special or coincidental cases that are possible) means in our 
case that all zeros of V* are simple and disjoint with A. 

For p = 5 (6 points) we can find essentially two different configurations: the linearly ordered set 
(Figure [s] left) and the branched tree (Figure [s] right). In general, a non-degenerate r(^) consists 
of a set of linear branches, like in Figure [7j left, and Figure [8j left, and it may have a branch point 
at a zero of V* where three branches merge, like in Figure |8j right. 

Thus, r*(^) is an analytic tree. It will be crucial for our study of families of positive ^-critical 
measures below. In fact, r*(^) plays a role of the "origin" of a coordinate system on the parameter 
plane of the above mentioned families of measures. 



Observe finally that by Lemma 5.11 Proposition 5.17 and the result of Lavrentiev, the Robin 
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Figure 6: Chebotarev's compact for three points. 

(or equilibrium) measure Ar of T = r(^) is a (unit) positive ^-critical measure with a connected 
support, and that Ar is the only (unit) positive ^-critical measure with this property. 

6.2 Cuts of minimal capacity and convergence of Pade approximants 

For A = {ao, . . • , flp} C C we denote by U{A) the class of analytic germs f{z) = Yln>o fnZ~^ at 
infinity that admit an analytic continuation to C\^. For a fixed germ / G U{A) consider the family 
5'(/) of cuts F G C which make / single- valued in their complement: 

5^(/) = {F CC: f holomorphic in C \ F}. 

We are interested in the set F = F(/) C d{f) such that 

cap(F) = min cap(F). (6.2) 

Observe that this problem is a generalization of that considered in the previous section. 

Approximately 40 years ago one of the hottest topics in approximation theory of analytic func- 
tions was the problem of convergence of diagonal Pade approximants [n/n]f = P/Q for functions 
/ G U[A). In this connection J. Nuttall made a basic conjecture that the sequence [n/n]/ converges 
(in capacity) to / G U{A) in C\F(/); he also proved his conjecture in some special cases (see review 
|,51j). The problem was completely solved in 1986 by H. Stahl, who proved Nuttall's conjecture in 
a striking generality: for closed sets A with cap(^) = 0. In particular, he established that for any 
/ G U{A) there exists an essentially unique set F(/) satisfying ( |6.2[); h e also characterized it in 



terms of the equilibrium conditions and the S'-property (see Section 5.3). We refer the interested 
reader to jMl CSl EH [Zl • 
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Figure 7: Chebotarev's compact for 4 points, when the zeros of V* are simple (left) or double. 




Figure 8: Chebotarev's compacts for 6 points. 

We note that our concept of an ^-critical measure can also be extended to compact sets A of 
capacity zero, but in this paper we deal with finite sets A only, assumption that we keep in the 
sequel. 

For a finite set A there obviously exists only a finite number of different possible solutions (sets of 



minimal capacity) of (6.2 ), that we denote by Fq, . . . , F^v, in such a way that Fq is the Chebotarev's 
continuum for A. In other words, for any / G h({A), F(/) € {Fq, . . . ,F7v}. 

For any k = 0, . . . , N , the equilibrium measure A = of F = F^ is ^-critical and satisfies 



U (z) = Pk = const, z gT = supp(A), (6-3) 
d d 

U^{z) = - — U^{z), z at regular points of F, (6.4) 



where n± are the normal vectors to supp(/i) at z pointing in the opposite directions. Moreover, it 



35 



follows from [70^ that these conditions define A and T uniquely in the given homotopic class (we 
spare details here). 

A comparison with Lemma 5.11| (for tp = 0) shows that among all ^-critical measures the Robin 
(or equilibrium) measures are distinguished by the equality of the equilibrium constants: if 
F/c 1, . . . , Ffc M are the connected components of F = F^, then 



U\ 



- k.l 



u\ 



Chebotarev's continuum Fq corresponds to M = 1; the other F^'s may be combinatorially charac- 
terized using the structure of Fq. We present an example explaining this statement. 

Example 6.1 Assume p = 3, uq = —x — iy (with x > y > 0), oi = Ziq, 02 = —clq, 03 = — oi, like 
in Figure [7j left. Then there are exactly two different compact sets of minimal capacity. One is the 
Chebotarev's set Fq (as the one depicted in Fig. [Tj left), corresponding for instance to 



{z - ao){z - ai)y/'^ 



{z - a2){z - as), 

and the other, Fi, with two components (as in Fig. [9]), corresponding to function 

'{z - ao)(z - oi) 



£U{A). 



{z - a2){z - as). 

Observe that Fi is topologically equivalent to Fq with the arc connecting and removed. The 



\ 



as 



Figure 9: A minimal capacity set for 4 points, 
quadratic differential associated with A = Ai has a double zero v at the center of the rectangle, so 

t- - 



that (5.24) takes the form 



w - U^{z) = Re 



dt. 



which is the Green function of the Riemann surface of the function = A{x) with logarithmic 
poles at ooi and 002. This is also an elliptic integral of the third kind. 
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For a general set A with an arbitrary number of points consider the Chebotarev's continuum 
To and select any arc /3 of To connecting two zeros of the corresponding quadratic differential (a 
"zero-zero connection"). Then there exists a compact Ti of minimal capacity with two connected 
components, which is topologically equivalent to To \ (3, whose quadratic differential exhibits a 
double zero instead of the pair of zeros we have selected. We can repeat this operation with any 
other remaining zero-zero connections, until all these connections are gone. If p = 2ni — 1, the 
maximal number of connected components of a set of minimal capacity is m. The unique T with m 
components is again the zero level of he Green function for the Riemann surface of = A{x). 

The compact sets and their Robin measures associated with A play a central role in any 
investigation of the strong asymptotics for complex orthogonal polynomials - denominators of the 
Fade approximant of functions from U{A), see [52]. 



6.3 Further connections 

The central part in the Stahl's solution of the convergence problem for Fade approximants was a 
new method of investigation (based directly on the 5-property) of the n-th root asymptotics of 
polynomials satisfying complex non-hermitian orthogonality conditions, like those verified by the 
denominators of the Fade approximants to functions from U(A) . This method was further developed 
in [24^ in relation with the best rational approximations; the n-th root asymptotics was obtained 
for complex orthogonal polynomials with respect to varying weights (i.e. depending on the degree of 
the polynomial) . The existence of a varying weight motivates the appearance of an external field in 
the associated equilibrium problem. Accordingly, the 5-property in the related existence problem 
should be modified to include the external field too (we omit here the non-essential details). 

Let (f he a harmonic function in a domain $7 C C. For a curve F C O let A = \r,ip be the 



unit equilibrium measure on F in the external field ip, so that conditions (2.8) hold. Recall that 
additionally F has the 5- property if at the regular points of supp(A), 

^U\z) = J^U\z), (6.5) 

where n± are the normal vectors to supp(/u) at z pointing in the opposite directions (it is assumed 
that the set of irregular points of supp(A) has capacity zero). We call the support of such a A an 
S-curve. The cornerstone to any application is the problem of existence of such a curve, that we 
discuss here briefly. 

Given a harmonic function in a domain 0, and a homotopic class ^ of curves F C ^, we should 
find a curve F G 5^ with the S-property. 

A direct and constructive approach to the solution of this problem is based on the observation 
that if such a curve F G 5^ exists, its equilibrium measure A = Xp^ip is {A, y5)-critical for some set A of 
fixed points, which depends on the definition of ^ and on the singularities of f. Then (see Theorem 
5.1 and Lemma 5.2) we conclude that {C^ + $')^ = R, where R is some function, meromorphic 



in 0, with (usually, known) poles at A and (usually, unknown) zeros. These zeros are the main 
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parameters of the problem; they must be found using a system of equations, typically in terms of 
periods of ^fR dt) , that reflect all the given information (including the geometry of the class ^) . 

For (p = this is basically a classical method, which goes back to Abel and Riemann (abelian 
integrals, see the discussion below). The existence of an external field does not change the nature of 
the problem, but it posses additional technical difficulties. If ^ 0, we generally have supp(A) T 



(compare with (6.3) and (6.4)), and finding the support of the equilibrium measure might turn out 
to be a formidable task even for a fixed T. We refer to [231 for further details. Se also |43], as well 
as [6], where the ^-problem for Q = C, and ip = Re(P), where P is a polynomial, was considered. 

Another way to prove the existence of the S-curve independently of its construction is based on 
the electrostatic interpretation of the critical measures, which yields the following extremal problem. 



Consider the equilibrium energy E^{-) (see (2.5 )-(pi3|)) of a curve F £ ^ as a functional on 



Under rather general assumptions it is possible to prove that if a curve F G 5^, satisfying 

E^[r]= max E^[F] (6.6) 



exists, then it has the S- property. For (f = Owe have E^[F] = — log (cap (F)), and (6.6) is equivalent 



to the minimal capacity problem considered above. For 99 ^ 0, this is the weighted capacity (2.7) 



minimization (see e.g. 161]), and the method was outlined in p3] in connection with the best rational 



approximation of exp(— x) on [0,+oo). The discrete analogue of problem (6.6) and its connection 
with Jacobi polynomials was discussed in |41] . 

In a rather surprising twist, a completely different problem was reduced to the existence of 
an S-cuTve in [33], where the semiclassical solution of the focusing nonlinear Schrodinger equation 
was constructed using methods of the inverse scattering theory. The problem itself, as well as the 
methods of its solution, did not have a priori any visible connection with those examined in |24] . 
A partial explanation of the mystery is suggested by the connection of the orthogonal polynomials 
with the inverse scattering via the matrix Riemann-Hilbert (RH) problem (see e.g. |15j and |18]). 
Apparently, the corresponding RH problems considered in [24, and [3J5j are similar. 

The recently developed tools of asymptotic analysis of the RH problems of a certain class, 
such as the non-linear steepest descent method of Deift and Zhou (see [131 El IE] i ^ well as ^ ) 
combined with the 9-problem ( |46| 14 7| ) . have become powerful weapons in the study of the strong 
asymptotics of polynomials of complex orthogonality. We can find multiple examples in a series of 
works of Aptekarev, Baik, Deift, Kuijlaars, McLaughlin, and Miller, to mention a few (see e.g. H], 
[15] . [35] . [36] . [37] . [38] : this is necessarily a very partial list). 

One of the main ingredients of the solution of such kind of asymptotic problems (independently 
of the approach we follow) is the analysis of the S-property related to the concrete situation. The 
total potential of the corresponding critical (or equilibrium) measure is called the ^-function (see 
[15]). This function typically accounts for the leading term of the asymptotics, and the support of 
the measure is the set where the essential oscillatory behavior takes place. Hence, the construction 
of the (7-function or of the S'-curve, or even establishing the existence of the latter without finding 
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the parameters explicitly, is an important problem. In the presence of a significant external field 
it has been solved so far for some particular situations. The use of the critical measures in this 
context may present a new approach to the problem at large. 



7 Weak limit of zeros of Heine- Stieltjes polynomials 

We return to our original motivation, armed now with the tools developed so far, in order to analyze 
the possible weak limits of the polynomial solutions of (1.3). We formulate first a statement slightly 
more general than necessary for our problem. 

We are given A = {ao, ai, . . . , ap} C C, and a sequence of external fields of the form 



Re^>n, -^niz) 



E 

k=0 



Pk{n) 



\og{z - ak) , 



(7.1) 



where Pk{n) £ C. 

Theorem 7.1 Let G SO^n, n G N, 6e a discrete {A, '^n)- critical measure corresponding to an 
external field (7.1). If for a subsequence C N, limits 



^^^r = Pk^ fc = 0,l,, 



,P, 



exist, then any weak-* limit point fi of the normalized measures {/x„/n}, n £ J\f , is a continuous 
{A, If) -critical measure with respect to the external field ip given by ( |5.1[ ). 
In particular, if 



1 



(7.2) 



A:=0 



then fi is a unit continuous {A, (p)- critical measure. 



Proof. Without loss of generality, we assume that Un = Pn/n -—>■ n G J\f, where —>■ means 
convergence in the weak-* sense. By Proposition 3.3, if (7.2) holds then supp(/i„), n G M, are 
uniformly bounded, so that the set of normalized measures f„ is weakly compact, and /i is a 
probability measure on C. 



By Remark 3.8 it is sufficient to show that for smooth functions h, condition 

d 



dt 



implies 



dt 



0, 



n G N, 



(7.3) 
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Let 



k=l 



''fc 



reasoning as in the proof of Lemma |3.7| we conclude that (7.3) is equivalent to the condition 

.2 



n 



_ An) 



(7.4) 



k=l 



Observe that 



E 



^(n) _ Jn) 



lim 



h{x) - h{y) 

\x—y\>t X ~ y 



dVn{x)dUn{y). 



Since — h{y))/{x — y) is continuous, standard arguments show that 



1 ^/.(cD-Mcr) 

hm ^ > 



/i(x) - /i(y) 



X - y 



dfi{x)dfi{y). 



On the other hand, ^Pn/n ^, n G M, locally uniformly in C \ „4, with ip given in (5.1 ). Since 
h vanishes on A, are continuous, and 

n^N n ^ — ' neA/ n J J 



In consequence. 



lim ^ 



E 



Jn) _ Jn) 



k=l 



with defined in (3.10). Using (7.4), we conclude that fip(n;h) = 0, and it remains to apply 
Lemma 13.71 □ 



We consider next a sequence of pairs (QnjVn) of Heine-Stieltjes polynomials Qn of degree n 
and their corresponding Van Vleck polynomials Vn- One of the central results of this paper is a 
description of all the possible limits of the normalized zero-counting measures of the Heine-Stieltjes 
polynomials Qn- Since the residues in (1.6) are independent of n, by applying Theorem 7.1 with 
if = we get: 

Corollary 7.2 Any weak-* limit point of the normalized zero counting measures ^{Qn)/^ of the 
Heine-Stieltjes polynomials is a unit continuous A-critical measure. 
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Remark 7.3 Taking into account Proposition 3.2 we could restate the last result in terms of the 



zero-counting measures of Heine-Stieltjes polynomials corresponding to a generalized Lame equation 



with coefficients depending on n. An analogue of Theorem 7.1 has been used in [35] for the study of 
the weak-* limits of the normalized zero counting measures of the Heine-Stieltjes polynomials with 
varying positive residues and ^ C M. 



By [66j, the zeros of Van Vleck polynomials accumulate on the convex hull of A, so that the 
set of all Van Vleck polynomials is bounded (say, in the component-wise metrics). Hence, in our 
consideration of a sequence of pairs {Qn, Vn) of Heine-Stieltjes polynomials Qn of degree n and their 
corresponding Van Vleck polynomials Vn we suppose without loss of generality that there exists a 
monic polynomial V of degree p — 1 such that 

lim Vn = V. (7.5) 

n^oo 



Theorem 7.4 Under assumption (7.5), the normalized zero counting measure v{Qn)/T^ converges 
(in a weak-* sense) to an A-critical measure /i G TWi; furthermore, the quadratic differential 

V . 



■w 



A 



[z) dz' 



is closed, the support T = supp(/x) consists of critical trajectories of w, C \ T is connected, and we 
can fix the single valued branch of jA there by lim^^oo z\/y{z) /A{z) = 1. With this convention, 



lim \Qniz] 



|l/n 



exp Re 




{t)dt 



(7.6) 



locally uniformly in C \ T, where a proper normalization of the integral in the right hand side is 
chosen, so that 

-{t) dt — log l^l 



lim I Re 

2— >00 




0. 



In other words, weak-* limits of the normalized zero counting measures of QnS are unit positive 
^-critical measures. The inverse inclusion (that any unit positive ^-critical measure is a weak-* 
limit of the normalized zero counting measures of Heine-Stieltjes polynomials) is also valid, but it 
cannot be established using methods of this paper. We plan to present the proof in a subsequent 
publication related to the strong asymptotics of Heine-Stieltjes polynomials. However, in the rest 
of the paper we identify both sets of measures. 



Proof. Let ^ be any weak-* accumulation point of {v{Qn)/'n}. By Lemma 5.2 the Cauchy trans- 



form of ^ must satisfy an equation of the form (5.4). Rewriting (1.3) in the Riccati form and taking 
limits as n — > cxD (with account of ( |7.5| )) we conclude that 

V 



{c^{z)y 



A 



z ^ supp(/i) 
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and r = supp(/i) is a union of critical trajectories of w. The rest of assertions about w follows from 
Proposition |5.17[ 



Finally, (7.6) is the straightforward consequence of (5.4) and the fact that for any monic poly- 



nomial P, log 1/|P(2;)| is the potential of its zero-counting measure. □ 



Theorem 7.4 provides an analytic description of the weak-* limits of the zero counting measures 
associated to the Heine-Stieltjes polynomials. However, this description is in a certain sense implicit, 
since it depends on the limit V of the Van Vleck polynomials Vn, that constitute therefore the main 
parameters of the problem. We must complement this description with the study of the set of all 



possible limits V . As it follows from Theorem 7.4 this can be done in two steps: 



i) describing the global structure of the trajectories of closed rational quadratic differentials with 
fixed denominators on the Riemann sphere, and the corresponding parameters (numerators). 
This problem has an independent interest; 

a) extracting from this set the subset giving rise to positive unit ^-critical measures. 

It was mentioned in Section |4] that problem i) is in general very difficult. In particular, we can find 



critical trajectories of any homotopic type. We have seen in Section 5.4 that not all of them will 
correspond to the support of an ^-critical measure, and elucidating this relation is the main step 
towards the complete description of the weak-* limits of the zero counting measures of Heine-Stieltjes 
polynomials. We start with a detailed discussion in the next Section of the simplest non-trivial case 
of three poles = 2), corresponding to Heun's differential equation. For a general p the geometry 
becomes so complex, that in this paper we just outline the main results (see Section [9]). 

8 Heun's differential equation {p = 2) 

In this section we concentrate on the differential equation ( |1.3[ ) with A{z) = [z — ao)(z — ai){z — 02) 
and V{z) = z — v. Observe that the Van Vleck polynomials constitute now a 1-parameter family, 
that makes the whole analysis much easier. So we introduce the quadratic differential 

V — z 2 

Wv = -^T—r dz (8.1) 

and two sets, 

V == {f G C : Wy is closed} , 

as well as the Van Vleck set 

V-f = {i) G C : ?; is an accumulation point of the zeros of Van Vleck polynomials} . 



A direct consequence of Theorem 7.4 is that 

V+ C V. 
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This inclusion is proper. Obviously, our main purpose, motivated by the analysis of the Heine- 
Stieltjes and Van Vleck polynomials, is to study V+; along this path related questions will be dealt 
with, such as the structure of the closed quadratic differentials of the trivial homotopic type, and 
the set of positive critical measures. On the other hand, general quadratic differentials and signed 
critical measures have an independent interest, and some results will be presented below. 



8.1 Global structure of trajectories 

Along this section we denote by Ty the critical graph of lUy, that is, the set of critical trajectories 
of together with their endpoints (critical points of zuy). 

There exists a unique v* = v*{A) E C such that the critical graph T^* is a connected set; 



Ty* coincides with the Chebotarev's compact T* associated with A (see Section 6.1). To simplify 
terminology, in the context of three poles we call v* the Chebotarev's center for A. 

Since the value v* is in many senses exceptional, along with the poles Oj, we will introduce the 
notation ^* = {v*} for the "exceptional set". 

We have mentioned in Section [4] that the global structure of the trajectories of a quadratic 



differential can be extremely complicated. For the quadratic differential (8.1) a certain order is 
imposed by the double pole at infinity with a negative residue. 

Proposition 8.1 Let A{z) = {z — aQ){z — ai){z — 02) he a polynomial with simple roots in C and 



V G C\.4*. Then the quadratic differential (8.1) has a closed critical trajectory (3 containing v. Let 



$7 he the hounded domain delimited hy /3. Then Q contains at least two points from A. 



Proof. Due to the local structure of the trajectories, has a closed trajectory freely homotopic 
to infinity (in other words, it is topologically identical to a circle and contains all the finite critical 
points of in the bounded component of its complement). According to Theorem 9.4 of [73], this 
closed trajectory is embedded in a uniquely determined maximal ring domain T^oo swept out by 
homotopic closed trajectories of ruy. We denote by j3* the bounded connected component of the 
boundary of T^oo- 

Obviously, (3* contains at least one critical point of tu^,, and no finite critical points of lie in 
7^00 (that is, in the unbounded component of C \ /?*). We conclude that all finite critical points of 
Wy lie either on (3* or in the bounded component of its complement. If this component is empty, it 
means that (3* contains all the critical points, and it is the Chebotarev compact for A. Otherwise, 
the bounded component of C \ /3* contains an interior point, and from the local structure of the 
trajectories of at simple poles we infer that [3* cannot contain only poles of vj^. Hence, v £ P* , 
and at least two of the three trajectory arcs emanating from v belong to /?*. It is easy to see that 
either they end up at respective poles (and then again (3* = T*), or they form a closed loop, that 
we call p. We call O the bounded domain delimited by /3. There is only one trajectory emanating 
from V remaining, that either is recurrent or ends at a pole from A. Hence, /3* cannot contain more 
than one pole, so that O contains at least two poles. This concludes the proof. □ 
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Remark 8.2 An examination of the proof of Proposition 8.1 shows that the existence of the ex- 



tremal trajectory /?* containing a zero of y such that all trajectories outside /?* are closed and 



homotopic to infinity is a fact, valid for an arbitrary quadratic differential of the form (4.1 ) 



Let V ^ C\A* . By the theorem on the local structure of the trajectories of zu^j, there are three 



trajectories originating at v. From Proposition 8.1 it follows that two of them form a single closed 



loop P that splits C into two domains: the bounded component of the complement, O, and the 
unbounded one, that we denote by T>. Let us denote by 7 the remaining trajectory emanating from 
V. Observe that (3 and 7 have a single common point, v; according to the relative position of 7 with 
respect to (3 we can establish the following basic classification for the critical graph of Wv- 

1) Exterior configuration: 7 \ {v} belongs to the unbounded component D of the complement of 



P (see Fig. 10 left). In this case v £V (the quadratic differential is closed). 

One of our main results, which we discuss later, is that v £ V+ if and only if F^ has an exterior 
configuration. 

2) Interior closed configuration: 'j\{v} belongs to the bounded component 0, of the complement 



of P, and 7 is finite (see Fig. uJn right or Fig. 11 right). In this case 7 is a critical trajectory 



joining v with one of the poles aj, and v £ V (the quadratic differential is closed). 
3) Interior recurrent configuration: 7 \ {v} belongs to the bounded component 0, of the comple- 



ment of /?, and 7 is not finite (see Fig. 11 left). In this case 7 (in fact, all trajectories in fl) 



is a recurrent trajectory, dense in 0,, and v ^ V. 

In order to prove that our classification exhausts all the possibilities it is convenient to single 
out the following simple statement: 

Proposition 8.3 The quadratic differential Wy is closed if and only if there exists one critical 
trajectory of vjy. 

Proof. Obviously, we only need to prove sufficiency. Assume that 7 is a critical trajectory joining 
for instance v and ai. In this case, taking into account the residue at infinity, we conclude that 




dt = 



A{t) 

if we integrate along a simple arc in 17 \ 7 joining both poles. This means that 




is a harmonic function on the Riemann surface TZ oi = (x — 02) (x — 03) with two cuts along the 
lift of 7 to TZ. Reasoning as in |24j we conclude that a zero level curve of this function connects 02 
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Figure 10: Critical graph of a quadratic differential in an exterior (left) and closed interior config- 
urations. 



and 03, and its projection on C constitutes the second critical trajectory of Since the critical 
graph r^, of is a compact set, the differential is closed. 

The remaining case is analyzed in a similar way, and this concludes the proof. □ 



Recall that by construction (see Proposition 8.1 1, /3 is part of the boundary /3* of the maximal 



ring domain T^oo swept out by homotopic closed trajectories of w^. Hence, in the case of an exterior 



configuration, 7 C /3* is a critical trajectory. By Proposition 8.3 w G V. An analogous conclusion 
is obtained if 7 \ \y\ C O is critical. 

Finally, assume that 7 \ \y\ C is recurrent. According to Corollaries (1) and (2) of Theorem 
11.2 in [71], its limit set is a domain bounded by the closure of a critical trajectory. Since in this 
case no critical trajectories can exist in we conclude that the closure of the limit set of 7 is Jl, 
which concludes the proof of the statements above. 

Let us summarize part of our conclusions in the following statement: 

Theorem 8.4 Let u G V \ .4*. Then the critical graph of is the union of three critical 
trajectories, 

= aU/?U7, 

such that (3 = f3{v) is a closed loop that contains v, a = a{v) is an arc in the hounded component 
i7 ofC\(3 joining two poles from A, and 7 = j{v) connects the remaining pole with the zero v. 
The set C \ F^, is the union of two disjoint domains: 

• the bounded component $7 \ F^ 0/ C \ F^ is a ring domain bounded either by a and (3 (for an 
exterior configuration) or by F^ (for an interior configuration); 
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Figure 11: Trajectories of a quadratic differential in an interior recurrent (left) and interior closed 
configurations. 



• the unbounded component 2? \ r„ = TZoo o/C \ r„ is a disc domain in C, bounded either by (3 
and 7 (for an exterior configuration) or by (3 (for an interior configuration) . 

For what follows we will fix an orientation of both arcs a and 7. For instance, we can agree that 
a goes from Oj to aj if i < j, and 7 always goes from a pole to v. 

8.2 Homotopic type of a closed differential 

The partition of the plane by the critical graph r„, associated to the closed quadratic differential Wv 



and described in Theorem 8.4 , allows us to introduce a geometric characterization of the trajectories 
of w "at large" . 

For E V we define the homotopic type of ro„ as the free homotopic class in C \ ^ of any of its 
closed trajectories in $7 \ F^. Any such a trajectory is a closed Jordan curve in C \ ^ containing 
exactly two points from A (say, ai and 02) in its interior, and leaving in the exterior the third pole, 
ag. We can further think that the (Caratheodory) boundary component given by the "two-sided 
a" belongs to this homotopy class (we think of the two-sided a as the closed curve with ai and 02 
in its interior). Then, without loss of geometric information, we can identify the homotopic class of 
the two-sided a with the homotopic class of a itself, considered now as an arc with fixed endpoints 
ai and 02 in C \ {oq}. 

Thus, for V & V\A* the homotopic type of Wv has a combinatorial component (namely, which 
two of the three poles from A are joined by a) and the geometric one (given by the homotopy of a 
in the punctured plane with the third pole removed) . 
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It follows from the general theory that there exists a one- (real) parametric family of closed 
differentials with the prescribed homotopic type. One standard way to parametrize such a 
family is by the tz7t,-lengths of the trajectories in Q \ Ft,; another one is based on the -aj^j-lengths of 
the conjugate (orthogonal) trajectories. More precisely, recall from Section |4] that the ti7^,-length of 
a curve r is 

„ „ 1 




By definition of the (horizontal) arc, given w S V, all closed trajectories of in \ have the 
same tu^j-length, equal to 

/ 9 1 1 ^ 1 1 

'"V — ^ 1 1 '^11 rot) ■ 

Moreover, is the minimum of the ro^-length of a closed Jordan curve separating the boundary 
components of O \ F^,; it is called the length of the circumferences of the cylinder associated with 
n\T^ (see IZi Chapter VI]). 
The conjugate value 

hy inf {llrllro^ : r connects boundary components of \ r^,} 

is called the height of the cylinder associated with $7 \ r„. Again, it coincides with the ro^-length 
of any arc of orthogonal trajectory connecting the boundary components of 17 \ 

Lemma 8.5 Let a be a Jordan arc lying in C\A (except for its endpoints) and connecting two 
poles from A. Then for any value h > there exists a unique v £ V such that Wy has the homotopic 
type a, and = h. 

This Lemma may be proved by reduction to a general Theorem 21.1 from |74| §21] on the existence of 
finite differentials. Our quadratic differentials are not finite (due to the double pole at infinity), 
but they can be approximated by finite differentials in a way preserving essential characteristics. 

The same family of differentials of the homotopic type a may be parametrized alternatively 
by the length 1^ of the circumferences of the cylinder associated with Q \ L^. In this case each 
homotopic type has a minimal (strictly positive) admissible length: 

Lemma 8.6 Let a he a Jordan arc lying in C\A (except for its endpoints) and connecting two 
poles from A. Define 

L = L{a) = mf{ly : u € V and has the homotopic type a]. 

Then L > and for any I > L there exists a unique v £ V such that zuy has the homotopic type a, 
and L, = I. 



Lemma 8.6 may be reduced to Lemma 8.5 or derived from the general existence theorems related 
to the moduli problem, see |711 §21]. 
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8.3 Correspondence between closed differentials and ^-critical measures 

We have proved in Section [5] that for any signed ^-critical measure fi with /i(C) = 1 there exists a 
closed quadratic differential in terms of which the measure and its potential may be analytically 
expressed. In general, this is not a bijection, since many critical measures correspond to the same 



quadratic differential. It follows from the proof of Proposition 5.17 that for p = 2 a one-to-one 
correspondence between closed differentials and signed critical measures is restored if we consider 
the ^-critical measures n with the additional property that C\supp(/i) is connected. This subclass 
is the most important for applications. 

Let V \ A*. Obviously, y^V/A, with V{z) = z — v, has a single- valued branch in C \ a U 7; 
the critical trajectories a = a{v) and 7 = j{v) were introduced in Theorem |8.4[ We fix the branch 
by requiring that lim^-^oo z^yV/A{z) = 1. Next we choose the positive (anti-clockwise) orientation 
in a neighborhood of infinity, that induces orientation on each side of a and 7. We denote by 
the subindex "-|-" the boundary value of a function at a and 7 from the side where the induced 



orientation matches the given orientation of each arc (see the remark after Theorem 8.4). 

With this convention, and taking into account that a and 7 are trajectories of Wy, we conclude 
that 

1 



defines a signed real measure on a U 7. Moreover, taking into account the residue at infinity, we see 
that J d^v = 1. Hence, we have proved the following 

Proposition 8.7 For any v £ V\^* there exits a unique signed A-critical measure fi^ with ^t,(C) = 
1, such that supp(/i„) = a U 7. Furthermore, fi^ is absolutely continuous on supp(//^) with respect 
to the arc-length measure, and formula (8.2) holds. 



Remark 8.8 This construction can be extended in a natural way to the Chebotarev's compact 
(v = V*) and to the degenerate cases when v £ A; in these situations, measure /ly* is positive. 
Hence, u 1— > //t, is a mapping from V into the set of signed unit measures on C. 

We introduce next an analytic function that will allow us to study the structure of the set V of 
points V £ C such that zuy is closed. 

Let vq £ V \ {v*}; recall that a = a{vo) joins the two poles of Wy^ not connected with vq by a 
critical trajectory. In a simply connected neighborhood of vq, disjoint with a, 




t - V 



is analytic in v, single- valued, and 

^'N) = i I dt + 0, (8.4) 

J a V - 'V^)A{t) 
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since this is a period of a holomorphic differential on the elhptic Riemann surface of the algebraic 
function = {t — VQ)A{t). This construction defines an analytic and multi-valued function w in C, 
with w' 7^ 0; however, formula (8.3) allows to specify a single- valued branch of it; in a neighborhood 
of a point only in C \ {v*}. 

Proposition 8.9 For every vq ^ V \ {f*} there exists a neighborhood B of vq such that V r\ B 
contains an analytic arc I passing through vq, and such that the homotopic class of the trajectories 
of for V € i is invariant. 

Proof. In a small and simply-connected neighborhood B if vq consider the branch of w given by 
formula (8.3), so that 'w{vo) = 7rz(l — fi{'y{vo))). Since w'{vq) / 0, the level curve 

£ = {v : Rew{v) = 0} 

is well defined in B, and constitutes an analytic arc passing through vq. Clearly, if f E BCiV is such 
that the homotopic class of the trajectories of Wy and the same, then necessarily v G i. 

Reciprocally, assume v £ i D B, and consider the trajectory of that emanates from the same 
pole as 7(^0)- Due to the continuity of the level curves of 




with respect to a variation of v, this trajectory is either critical (and then the proposition is proved) 
or recurrent. In the latter case it must intersect the orthogonal trajectory of Wy starting from v 
(see |74| §11]), which contradicts the hypothesis that Yievu{v) = 0. □ 

Now we can describe completely the structure of the set V: 

Theorem 8.10 The set V is a union of a countable number of analytic arcs ik, k £ Z, each 
connecting v* and 00. 

Two arcs from V are either identical or have v* as the only finite common point. The homotopic 
type of the critical trajectories of in C\A remains invariant on each arc Ik \ ■A.* . 
There are three distinguished arcs i^j k £ {0,1,2}, such that 

(i) Ik connects v* with infinity and passes through a^; 

(a) for every v G £k the homotopic class of trajectories of the closed quadratic differential Wy is 
trivial; 

(Hi) Function ^^{^{v)) is monotonically decreasing from iJ-v'iliv*)) to —00 as v travels Ik from v* 
to 00. 
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Proof. Assume that Oj's are not collinear (the analysis of the colhnear situation is simpler). For 
vo = ao € V, the arc 7(^0) vanishes, and q(?;o) = [01,02] is the straight segment joining ai and 02, 
so that we can fix the single-valued branch of tt; in a neighborhood of oq by 

w{ao) = 1 ^ — dt = iri. (8-5) 

Denote by 7i the half plane containing oq and determined by the straight line passing through 
oi and a2- Then 



w{v) 



ai,a2 




and (8.5) determines the single- valued branch of w in 7i. Let Tq be the level curve 

To = {z G Rew{z) =0}. 



We have that 

W{V*) = TTl (r* \ 7(1;*)) 

(r* \ j{v*) is the union of two arcs of the Chebotarev compact joining v* with oi and 02). Hence, 
V* G To. 

Since the rotations and translations of the plane do not affect the character of the level curves 
of w, we can assume that both ai, 02 S M. Then it is immediate to see that Tq can intersect M at a 
single point (which belongs to the segment [oi , 02] ) • The other end of Fq must diverge to infinity. 
This establishes the existence and properties of the distinguished arcs ik, k £ {0,1,2}, described 
above. 

Consider now the analytic function w in the infinite sector delimited by two contiguous distin- 
guished arcs £ki k G {0, 1, 2}. Fix there a vq £ V and take the single-valued branch determined by 
the condition 



w{vo) = J ^ ^y^^^j dt = Trifj,^^{a{vo)). 

Then the level curve i = {Kew{v) = 0} is an analytic curve passing through vq, that can intersect 
the boundary of the sector only at v* . Hence, i joins v* with 00. The number of different curves I 
is given by the number of different homotopic types of closed trajectories, which is countable. This 
concludes the proof. □ 

Remark 8.11 As v approaches v* along an arc ik C V, the support supp(/x„) = aUj tends to 
the Chebotarev set F*, but possibly covered several times, in accordance with the homotopy class 
of vjv on ik- 

Finally, it is convenient to consider another independent parametrization of measures fiv, v G V, 
in order to connect the characteristics of their logarithmic potential U^^ with the geometrically 
defined values of the corresponding quadratic differential zu^. Applying Lemma 5.11 we get 
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Lemma 8.12 Any measure fi^, v G V, is characterized by the following property: U^^^ is constant 
on each connected component o/supp(^), 

[/^"(z) = Ca for z G a, U^^''{z) = for z G 7, 

and at any regular point o/supp(/i^), 

^U^'-[z) dU^'-^z) 



dri- 



^.6) 



where n± are the normal vectors to supp(;U^) pointing in the opposite directions. Moreover, we have 
the following relations: 



l^ = 2fi^{a) = 2 (1 - Hvh)) = 2iriw{v), K = A^a 
8.4 Positive ^-critical measures 



7r(cQ, - c^). 



Theorem 8.10 gives a complete description of the set V of points v that make the quadratic differen- 
tial in (8.1) closed. By Proposition 8.7, to every v gV \t corresponds a unique signed ^-critical 
measure given by formula (8.2). Our next goal is to isolate the subset 

V+ = {w G V : /ii, is positive}. (8-7) 

Proposition 8.13 Let w G V. Measure fi^ is positive if and only if either v = v* or v is in an 
exterior configuration. 



See Section 8.1 for the definition of the exterior configuration. 

Proof. The case v = v* is trivial, so let us assume that v ^ v* . Consider 



u{z) = Re 



It-v 



dt , u{z) = Im 



It-v 



dt. 



Since a = a{v) and 7 = 7(f) are trajectories of lu^, function u is single- valued and harmonic in 
C \ r, continuous up to the boundary, and the closed trajectory (3 (see Proposition 8.1 ) is its zero 
level curve. By the selection of the branch of the square root, u{z) ~ log \z\ as z — > 00, and we see 
that u{z) > for z G C \ ($7 U 7), and u{z) < for 2; G O \ P. In consequence. 



dn 



u{z) > on a. 



and on 7, 



dn 



u{z) 



> 0, if u is in an exterior configuration, 
< 0, if w is in an interior configuration. 
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where d/dn denotes the derivative in the sense of the outer normals. By the Cauchy-Riemann 
equations, 

-7ru[z) = -7z-u[z), 
OS on 

where d/ds is the derivative along each shore of the cuts a and 7 in the direction of the induced 
orientation. Hence, we conclude that /^?;|^ is always positive, while ^iv\^ is negative if and only if v 
is in an exterior configuration. □ 

Remark 8.14 Observe that we have proved that always piv\^ > 0. For v v*, hy construction 
fiv{a) + fJ-vil) = 1, and does not change sign on each connected component of T, so that fiy is 
positive if and only if /i^(a) < 1. An equivalent condition can be stated in terms of the ro„-length 
of the critical trajectories a and 7: 

fJ'v>0 ^ ||a|U„ + ||7lU„ = 1- (8.8) 

Remark 8.15 Figure[3]illustrates that only in the exterior configuration the ro^-rect angles intersect 
the support of fi^ only once (cf. Lemma 5.3). 



Proposition 8.13 provides an "implicit" geometric description of the set V+. Our main result of 



this section describes this set completely: 



Theorem 8.16 Let Ij^, k E {0, 1, 2}, be the distinguished arcs in V described in Theorem 8.10 The 
set V+ is the union of the sub-arcs £^ of each £k, k G {0, 1, 2}, connecting Uk with the Chebotarev's 
center v* (and lying in the convex hull of A). 

Furthermore, let us denote ruk = ^t,*(r^), k G {0,1,2}, where is the arc of the Chebotarev 
compact T* connecting v* with (mo + mi + m2 = Ij- If v £ ik <^ V+, k £ {0, 1,2}, then the 
trajectory 7(f) connects v with the pole ak and 

< l^viliv)) < ruk. (8.9) 

In this case both trajectories j^v) and a{v) are homotopic to a segment. The bijection fJ-viliv)) v 
is a parametrization of the set n V+ by points of the interval [0,mfc]. 

Proof. Straightforward estimates show that v fiy{a{v)) is unbounded on each arc ik C V. 
Furthermore, 

IJ,v{a{v)) = 1 ^ ^^{'j{v)) = ^ V e A. 



Since always fiy{a{v)) > (see Remark 8.14), we conclude that iiy{a{v)) takes values in (0, 1) only 
on the portions of the distinguished arcs £k, k G {0,1,2}, joining the Chebotarev center v* with 
each pole. For any other arc ik, ^„(a(t;)) > 1, and fiy is not positive. □ 

See illustration of the correspondence between the position of w on V and the trajectories of ro„ 
in Figure [12 
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Corollary 8.17 Any positive A-critical measure /x„ has the trivial homotopy type. 



Remark 8.18 Taking into account this Corollary and the definition of we can rewrite (8.9) is 
the following equivalent form: 



where j'l = min({0, 1, 2}\{/c}), j2 = max({0, 1, 2}\{fc}), and we integrate along the straight segment 
joining aj^ and • This system of equation defines the set V+ completely. 

Finally, in relation with our primary goal of the description of the weak-* asymptotics of the 
zeros of Heine-Stieltjes and Van Vleck polynomials we state the following important result, which 
however will not be proved completely in this paper. 

Theorem 8.19 V+ = V+. 



The inclusion V+ C V+ follows from the definition of V+ and Theorem 7.4 if f G V+, then fi^ 
is a limit distribution of the zero counting measures of the Heine-Stieltjes polynomials, so that 
fiy > 0. The inverse inclusion {fi^ > =^ v € V+) is also valid, but it cannot be established using 
methods of this paper. We plan to present the proof in a subsequent publication related to the 



strong asymptotics of Heine-Stieltjes polynomials. However, as the consequence of Theorem 8.19 



we identify the set of accumulation points of the zeros of Van Vleck polynomials V+ with the values 
of V £ V making fj,^ > 0. 

Remark 8.20 Although V+ and the Chebotarev compact T* are topologically identical and, ac- 
cording to numerical experiments carried out by B. Shapiro, metrically very close, they are not the 
same (an in consequence, the conjecture made in [67J is false). For simplicity, take ao = 0, oi = 1, 
Ima2 > 0, and define for z in the upper half plane close to the origin 





/i(^)= / \h^dt, h{z) 



where we integrate along segments joining with z. The change of variables t^zu'm the integrand 
of /i yields the asymptotic expansion 



On the other hand, using the asymptotic expansion of the integrand of /2 we get 
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where we take the main branch of \[z. 

Observe that Im /i = defines locally the set V, while Im /2 = corresponds to the Chebotarev 
compact of A. Assuming that both curves are tangent at the origin we conclude that v* j ^Ja^ > 0, 
so that V* lies on the bisector of the interior angle formed by 1 and 02 at the origin. In order to check 
the second order tangency, we can invert the mapping y = fi{z) and analyze F{y) = f2{fi^{y)) 
that maps the real line into itself at the origin. Setting v* = s(l + 02), s > 0, we get 

Setting V* = 8^/02, s > 0, we get 

(02 + 1)^;* f ^ 1 \ 
— = [^+ — ] s, 

which is real only if \a2\ = 1. This shows that at oq both curves V and T*, however close, are not 
identical, at least when the triangle with vertices at A is not isosceles. 



9 General families of ^-critical measures 

Most of the arguments presented in Section |8] for the case p = 2 may be carried over to the case 
of an arbitrary p with minor modifications. However, in a certain sense the multidimensional case 
is significantly more complicated. The volume of this paper does not allow to develop the whole 
theory, covering both signed ^-critical measures and closed quadratic differentials of an arbitrary 
homotopic type. In turn, without such a theory it is more complicated to separate positive ^-critical 
measures from the signed ones. So, we restrict ourselves here to a less ambitious goal allowing a 
shorter treatment: we put forward a constructive characterization of the positive ^-critical measures 
by proving that the constructed measures are indeed positive. The complete proof of the fact that 
there is no other positive ^-critical measures is matter of a forthcoming paper. 

9.1 Mappings generated by periods of a rational quadratic differentials 

Let us recall the notation. We have the fixed set A = {ao, ai, • • • , Op} of distinct points on C, 

V{z)'^ll{z-v,), R{z)'^^ and w{z) = -R{z) (dz)^ (9.1) 

is a rational quadratic differential on the Riemann sphere C. Zeros fj's of V are not necessarily 
simple, and we denote v = {vi, . . . , fp-i} G C^^-*^ with account of their multiplicity. Let also 

V = {v : w is closed}. 
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Occasionally, it is convenient to consider as a differential form ^\/li{z)dz on the Riemann surface 
TZ of or equivalently, on the hyperelliptic surface of genus p — 1 given by w"^ = A{z)V{z). 

Let r = 71 U • • • U 7p be a set consisting of p disjoint arcs 7^, each one connecting a pair of 
points from ALI v in such a way that C \ T is connected and ^/R is holomorphic in C \ F. The 
Caratheodory boundary of C \ F consists of p components % == 7^ U 7^ , with a positive orientation 
with respect to C \ T. We can consider % as cycles in C \ F enclosing the endpoints of 7^. Part of 
TZ over C \ T splits into two disjoint sheets, so we may consider % as cycles on TZ. 

Let us define 

Wk{v) = Wk{v,T) = — (j) y/R{z)dz, k = l,...,p, (9.2) 



Ik 



where V^];:^^ are the boundary values of the branch of \/^ in C \ r defined by lim^^oo R{z) = 1. 
Clearly, the boundary values {VR)± on 7^ are opposite in sign. Therefore, with any choice of 
orientation of 7^ and a proper choice of = (\/i2)+ on 7/5, we will have 



Wk{v) = — [ y^R{z)dz, k = l,...,p. (9.3) 

By the Cauchy residue theorem we have that wi + ■ ■ ■ + Wp = 1 for any v S C^^"*^. Thus, we can 
restrict the mapping v t-^ w to p — 1 components of w = {wi, . . . , Wp-i) G C^^-*^. In this way, we 
have defined the mapping 

V{-,T) : CP^i ^ CP-^ such that V{v,T) = w{v,T). (9.4) 

Each component function Wj{vi, . . . ,Vp-i) is analytic in each coordinate (even if is at one of 
the endpoints of 7^). Once defined by the integral in (9.3), this analytic germ allows an analytic 
continuation along any curve in C \ .4. Arcs 7^ are not an obstacle for the continuation since the 
integral in (9.3) depends only on the homotopic class of T in C \ (.4U i»). The homotopy of T is a 
continuous modification of all components simultaneously in such a way that they remain disjoint in 
all intermediate positions. Under this assumption we can continuously modify the selected branch 
of in C \ r along with the motion of T. 

We note that this notion of homotopy is different from the concept of homotopic class based 
on a choice of a collection of Jordan contours in C \ „4, which is standardly used to classify closed 
differentials (see |i74] ) . 

The homology basis {7j}jZi[ of the domain C \ T C 7^ defined above may be completed (in 

several ways) to form a homology basis 71, . . .7p-i, 61, ... , (5p_i for TZ. We can select the cycles 6j 
as a lifting to 7^ of a collection of arcs 6j C C \ T, each connecting new (different) pairs of points 
from AUv (see Figure 13 ). We denote by A = {61, ... , (5p_i), A = {5i, . . . , 5p^i). Accordingly, we 
define mappings 

Wk{v) =Wk{v,A) = ^ (f ^JR{z) dz, k = l,...,p-l, (9.5) 
V{;A) -.C^-'^ ^CP-^ such that V{v, A) = {wi{v, A), . . .Wp-i{v, A)). (9.6) 
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a2 



72 



71 



Figure 13: An homology basis for TZ. 

An important new mapping associated with the complete basis of homology on IZ is 

<P(-,r,A) : C^-i — >m2p-2 

given by 

<P(iJ,r, A) = (Imw;i, . . . ,Imzi;p_i,Im5;i, . . . ,Im{t;p_i) = (Im(P(r>, T)), Im(P(i>, A))). (9.7) 

In order to prove that both mappings, V and *P, are locally invertible we need the following lemma, 
which is a standard fact of the theory of the Riemann surfaces, see e.g. |68| Chapter 10] or |32^ 
Chapter 5]. 

Lemma 9.1 Let Wi, . . . , Wp-i be a basis of holomorphic differentials onlZ (a cohomology basis). 
Then 

(i) p — I vectors 



are linearly independent over C. 
(a) the system of 2p — 2 vectors, which consists of p — 1 vectors in (i) above plus p — I vectors 



D-l 



, j = l,...,p-l, 



are linearly independent over M. 

Proposition 9.2 Both mappings w = V{v) and {lm.w,lm.w) = ^{v) are locally invertible at any 
v{vi, ... , Vp-i) G (C \ A)P~'^ with Vi / Vj for i / j. 
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Proof. We have that for any j, k £ {1, . . . ,p — 1}, 



where 



dvk 



Vkit) 



1 

v{t) 

t-Vk 



dt 



1 



^/Rit) ^ . ^ , 

7, t-Vk 47r« J^^ y^A{t)V{t) 



Vk{t) 



dt 



Wk, 



and Wk 



1 



Vk{z) 



: dz, A; = 1, 



7j 



27ri ^A{z)V{z) 

Since Vfc are hnearly independent polynomials of degree p — 2 {V has simple roots), is a 

(i), vectors {dw/dvk}l.Z\ are independent; 



).l 



basis of holomorphic differentials on TZ. By Lemma 
this me ans that the Jacobian of V does not vanish and V is invertible. In a similar fashion, by 
(a), vectors {dw / dvk}^Z^\ and {dw/dvk}^Z^\, k = 1, . . . ,p — 1, are linearly independent 



Lemma 
over M. 



9.1 



The matrix of the linear mapping 
{Rev,lmv} G I 



p2p-2 



p2p-2 



is therefore nonsingular, and ^ is locally invertible. 



□ 



Proposition 9.3 Let fjP be an A-critical measure such that T = supp(/i^) has connected components 
7j^, . . . ,7p, and C\r is connected. Let vjq = RQ(z)(dz)^ be the quadratic differential associated with 
/u", where Rq = Vq/A, and = {vi, . . . ,Vp_^) is the vector of zeros of Vq. Assume that wq is in a 
general position (that is, all 's are pairwise distinct and disjoint with A). Then for an e > and 
any mj £ M., j £ {1, . . . ,p — 1}, satisfying 



Irrin 



1, 



there exists a unique solution v £ V of the system 



Wj{v,T) 



m 



l,...,p-L 



The quadratic differential w = —R{z){dz)'^ , R{z) = nfc=i(-^ ~ Vk)/A{z), is closed. The associated 
A-critical measure /i with supp(/i) homotopic to supp(/io) satisfies fi{-jj) = m 
where supp(/i) = 71 U • • • U 7p and supp(/x) is homotopic to supp(/i°). 



(9.8) 



■J' 



,P 



This result is a multidimensional version of Proposition 8.9 



Proof. We use the continuity of the dependence of short or critical trajectories of zu from the zeros 
-i) of V. Let, for instance, 7^^ be an arc of supp(//°) connecting two ^''-points Vi 



and For e > small enough the solution of (9.8) has two points, vi and f2, close to and ^2, 
respectively; moreover, one of the trajectories of w comes out of vi in the direction close the the 
direction of 7*^ at Wj*. This trajectory will be close to 7 and will pass near V2. If it does not hit vq 



then ^ f^^ ^/Rdt ^ M, in contradiction with (9.8). Same (even simpler) arguments apply to any 
trajectories connecting two points from ^ or a point from A with another from V. □ 
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This result introduces a topology on the set of ^-critical measures in a general position. We will 
call a cell any connected component of this topological space. By Proposition |9.3[ a cell is a manifold 
of real dimension p — 1. We can use as local coordinates either the vector (^ui = /^(7i), • . . , = 
Klp-i)) or V = {vi,.. .,Vp-i). 

Example 9.4 Consider the case p = 2, studied in detail in Section |8] In this situation the real 
dimension of each cell is 1. The image of each non-distinguished cell in the v-plane is an analytic 
are connecting the Chebotarev's center v* and oo. The unit positive measures, parametrized by 



the set V+ (see (8.7)), are represented on the f -plane as a union of three cells - arcs ij , j = 0,1,2; 



their boundaries are points from A* = {ao, ai, a2,v}, see Theorem 8.16 



For an arbitrary p, the boundary of a cell consists of pieces of manifolds of dimension < p. We 
come to a boundary point of a cell if either: 

(i) one of the components - arcs 7j C supp(/i), - degenerates to a point; this case is in turn 
subdivided into two subcases: 

(a) coalescence of a zero and a pole joined by an arc; reduction of p; 

(b) coalescence of two zeros joined by an arc; the polynomial V gains a double zero. 

(ii) two components (arcs) of supp()u) meet (at a zero of V). 

We need to take a closer look at the case (ii). Aga in, a simple but important example is p = 2, 
Section [sj An arc if G V+ is defined by (see Theorem 8.16) 



wi{v) = — / y/R{t) dt = £ (0,mi) C 



with a selection of the proper branch of the function. This is a "P-parametrization of the cell if , 
which uses the coordinates /ii = 'w{v,T), T = [ai,v]. The extremal values /ii = and fii = mi 
represent the boundary of the cell. Function w{v) is analytic at both points; it is convenient to 
analyze the reconstruction of fiy near = 0, that is, v = ai. The boundary point /ui = mi, 
corresponding to v = v* , is better seen from the point of view of the *P-mapping. 

Let 6 be an arc from 02 to v (homotopic to the arc [02, f*] C F* for v close to v*). Together 
with w{v) = w{v,'y) we consider function 

w{v) = — [ jR{t)dt. 

Then v = v* \s uniquely defined by equations 

lm.w{v) = lm.w{v) = 0, 

and moreover, a nearby point v on the arc lm.w{v) = is uniquely determined by a coordinate 
h{v) = Imt(;(t>) G M. It is important that h{v*) = 0, and h{v) changes sign when v crosses v* along 



the arc (see again Figure 12 where the reconstruction of supp(/i„) in dependence of v is illustrated). 
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At this moment it is not really important to determine which sign of h{v) corresponds to £^ . 
It is more convenient to introduce a number s = ±1 (depending on the branches of w and w) such 
that 

£^ = G C : lmw{v) = 0, lmw{v) = sh, h > 0}. 

This parametrization introduced originally around v = v* , may be then extended to the whole arc 
if. We note that h is the "height of the cylinder" (see [73]) and £ = 2fi = Ilew{v) is the "length of 
the circumferences", discussed in more details in Section 18.21 



9.2 Structure of the set of positive ^-critical measures 



6.1 



Let r* = r*(^) be the Chebotarev's continuum associated with A = {oq, . . . ,ap}, see Section 
it consists of critical trajectories of = —R*dz^, R* = V* /A. We will assume again a general 
position for the set A, that is, A and V* do not have common zeros and V*{z) = 111=1 (-^ ~ ''^k) 
does not have multiple zeros. Thus, the critical set A* = AU {v^, . . . , v*_-^} consists of 2p different 
points, and T* is comprised of 2p — 1 arcs, that are critical trajectories of zu. Each trajectory joins 
two different points from A*. 

We begin the c onst ruction of positive critical measures by introducing local i?-coordinates. As 



above (see Section 9.1), we identify measures /U G V+ with the corresponding polynomials V{z) = 
Y[^=\{z — Vk); furthermore, we define V by the vector v = {vi, . . . ,Vp-i) G C^^-*^ of its zeros 
(numeration is not important). Then each cell in V+ is a subspace of C^^"*^ ~ M^^^^, which is a 
manifold of the real dimension p — 1, defined by p — 1 real equations of the form 

lmwj{v)=lm(^-^. ^R{t)d?j =0, j = l,...,p-l, (9.9) 

where T = 71 U • • • U 7p-i U 7p is a union of Jordan contours % on TZ (double arcs) depending on v, 
but mutually homotopically equivalent for values of v from the same cell. Practically, any vq £ G 



has a neighborhood of v satisfying (9.9) with constant 7fc's. 

Now we come to the procedure of selection of combinatorial (rather than homotopic) types 
of cells; once the combinatorial type is fixed, the homotopic one will be determined from the 
Chebotarev's continuum, as described next. We start with the Chebotarev's continuum T* and the 
corresponding polynomial V*{z) = Y[k=i{^ ~ ^fc)- Each zero v* = v^, /c = 1, . . . ,p — 1, is connected 
by component arcs of T* with three other points, say a^, 02, 03 G A* . We select one of these three 
arcs (for definiteness, [u*,a|]) and join two other arcs to make a single arc [03,03], bypassing v* 
(we think that the arc [02)^3] still follows the two arcs from T* , but without touching v* , instead 
passing infinitely close to it). This procedure, carried out at each zeros of V*, creates a compact 
set r, and consequently, a cell G{T) of corresponding measures // G V+. 

The selection of T, and hence, of the cell G(r), is made by choosing one of the three connections 
for each there are 3^"^ ways to make the choice. Any choice splits T* into p "disjoint" arcs 
r* = 71 U • • • U 7p; out of them we select p — 1 arcs to make an homology basis for C \ T*, say 



71, . . . ,7p-i, and then consider the corresponding cycles as described in Section 9.1 
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Next, together with the hnk 7^ connecting to one of its neighbors from A* we will mark one 
more arc 6k connecting v^. with a different neighbor (i.e., a different point from A* connected by 



a branch to v^, see Section 



6.1). The choice of 7^ was arbitrary for each k; the choice of 
has to be made in such a way that p — 1 arcs jk and p — 1 arcs Sk are all different. Then the 
corresponding cycles {T, A} = {71, . . . 5i, . . . , 5p-i} form an homology basis on TZ and may 
serve to define the mappings V and ^ as in Section |9.1[ 

We will mention first the description of the cell G(T) in terms of the mapping V. This way of 
parametrization is equivalent (this equivalence is, however, not completely on the surface) to the 
"length of the circumferences" parametrization of the closed differentials (see Section |8] for the case 
p = 2). Let w = w{v) = V{v,T). We claim that the cell G(r) is completely defined by the system 

Wj{v)=tjeR+, j = l,...,p-l; (9.10) 

more precisely, there exists a domain M(r) = {(ti, . . . , tp_i) G ^} such that for any point 
(ti, . . . ,tp-i) G M(r) system (9.10) has a unique solution v £ C^~^. Moreover, the corresponding 



measure fi = fi^ satisfies /i(7j) = tj, and supp(/x) = 71 U • • • U 7^ = Ti, is homotopic to T. 

Summarizing, a rough description of the set V+ of unit positive ^-critical measures may be 
made as follows. The set V+ is a union of 3^~^ of closed bounded cells G{T) (T = 71 U • • • U 7p may 
be selected in 3^~^ ways). The interior G{T) of each cell consists of measures /i in general position 
with supp(/i) homotopic to T. Interiors of different cells are disjoint. Chebotarev's measure /u* 
(Robin measure of T*) is the only common point of all boundaries: /x* = Hr dG(T). 

The detailed proof of the assertions above and further analysis is beyond the scope of this paper. 
We will prove only that there exists a cell G(r) with the homotopic type T consisting of positive A- 
critical measures. It is easier to do it using the V- mapping (equivalent to the "height of cylinders" 
parametrization of the closed differentials). 

We consider the mapping {Im w, Im w} = ^{v; T, A), described in Section [9.1[ which is invertible 
in a neighborhood of v*. We select a vector (si, . . . , Sp-i) of signs: each Sj G {—1, +1}, and consider 

lmwj{v)=0, lmwj{v) = Sjhj, /i^ G M, j = l,...,p — 1. (9.11) 



For hj = 0, system (9.11 ) has a unique solution v* = (v^, . . . , v*^^). For sufficiently small hj > 
this system is still uniquely solvable. Equations Im Wj{v) = imply that differential zu in (9.1) is 



closed, the associated measure ^ is ^-critical, and supp(/x) = F^ = ji^^ U • • • U jp^v 

The homotopic type and signs of the components of fj, depend on the behavior of trajectories of 
w, which are originated at the points a* C A* and close to trajectories 6j. Any such a trajectory 
will hit the corresponding point vj if hj =0. If hj > 0, then it will pass from the left of Vj or from 



the right of Vj, see Figure 14 



A change from Sj to —Sj will change the direction of the turn. Therefore, there is a unique 
selection of vectors (si, . . . , Sp_i) such that all turns are right. Then the branch of '\/~R in C \ F^ 
will be close to the branch of in C \ F*, and therefore the corresponding measure /i will be 
positive. In this sense, the cell we entered contains some positive measures. Therefore, they are 
positive, since supp(^) are all homotopic. 
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Figure 14: Left and right turns. 
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